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Abstract. We describe an efficient algorithm to compute a pseudotrian- 
gulation of a finite planar family of pairwise disjoint convex bodies pre- 
sented by its chirotope. The design of the algorithm relies on a deepening 
of the theory of visibility complexes and on the extension of that theory 
to the setting of branched coverings. The problem of computing a pseu- 
dotriangulation that contains a given set of bitangent line segments is also 
examined. 
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1. Introduction 

1.1. Main result of the paper. Throughout the paper we address the prob- 
lem of computing efficiently a pseudotriangulation of a finite planar family of 
pairwise disjoint convex bodies presented by its chirotope : Here the term chi- 
rotope refers to a natural extension to finite planar families of pairwise disjoint 
convex bodies of the classical notion of chirotope (or order type) of a finite 
planar family of points OE]; and the term planar refers to any oriented topo- 
logical plane on M^, e.g., Euclidean plane, hyperbolic plane, Moulton planes, 
arc planes, etc.; cf. Appendix 1X1 

1.1.1. Chirotopes. Recall that the chirotope of a finite planar family of points 
is (or can be defined as) the map that assigns to each ordered triple of distinct 
indices of the family of points the position vector of the corresponding ordered 
triple of points, that is, the boolean vector of truth-values of the five relations 
"the third point of the triple belongs to the open left side (open right side, 
initial part, median part, final part) of the directed line joining the first point 
of the triple to the second point of the triple." Figure [T] shows five families of 
three points realizing the five possible chirotopes on the indexing set {1, 2, 3}. 
In this figure the plane is represented by the interior of a circular diagram. 
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Figure 1. 

marked with a little oriented circle to indicate its orientation, and each dia- 
gram is labeled at its left bottom corner with a symbol to name it and at its 
right bottom corner with the position vector of the ordered triple of points 
corresponding to the ordered triple of indices 1,2,3. The notion of chirotope 
of a planar family of pairwise disjoint convex bodies is defined similarly: as for 
families of points we use the notion of position vector as a coding of the relative 
positions of the convex bodies with respect to a line. To set out the definition 
we use the following standard terminology: a directed bitangent joining an 
ordered pair of disjoint convex bodies is, as illustrated in the left part of Fig- 
ure [21 classified left-left, right-right, left-right or right-left depending on which 
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sides (left or right side) of the bitangent are the convex bodies; walking along 
a directed bitangent we traverse successively, as illustrated in the middle part 
of Figure [21 its initial, median and final parts; the median part of a bitangent 
is called a bitangent line segment thereafter. Using this terminology we are 




Figure 2. 

able to define the chirotope of a finite planar family of pairwise disjoint convex 
bodies as the map that assigns to each ordered triple of distinct indices of 
the family of bodies the position vector of the corresponding ordered triple of 
bodies, that is, the boolean vector of truth-values of the twenty relations "the 
third body of the triple intersects the open left side (open right side, initial 
part, median part, final part) of the left-left (left-right, right-left, right-right) 
directed bitangent joining the first body of the triple to the second body of 
the triple." For example, consider the family of three convex bodies on the 
indexing set {1,2,3} depicted together with its 3 x 4 bitangents in the right 
part of Figure [2] (6 of the 12 bitangents are tritangents). Then its chirotope is 
the map x defined by 

right-left right-right left-left left-right 
X(l,3,2) =x(3, 1,2) = 11010 01010 10010 11010 
X(l,2,3) =x(3,2, 1) = 10000 01001 10001 01000 
X(2,3, 1) =x(2, 1,3) = 01000 01100 10100 10000. 

The number of chirotopes of planar families of 3 pairwise disjoint convex bodies 
on a given indexing set of size 3 is 531 and among these 531 chirotopes 118 
are simple chirotopes, that is, chirotopes of families of convex bodies with no 
tritangent; as for the chirotope of a planar family of points a key feature of 
the chirotope of a planar family of pairwise disjoint convex bodies is that it 
encodes its dual arrangement, i.e., the arrangement, in the space of lines of the 
plane, of the curves of tangents to the bodies; cf. Appendix El Throughout 
the paper we will assume that the boundaries of the bodies are free of line 
segments and that there is exactly one tangent through each boundary point; 
these assumptions facilitate the geometric definition of pseudotriangulations 
without ruling out any chirotope of families of pairwise disjoint convex bodies. 

1.1.2. Pseudotriangulations. Let 01,02,..., o„ be a finite planar family of n 
pairwise disjoint convex bodies; a boundary bitangent line segment is a bitan- 
gent line segment of the Oj contained in the boundary of their convex hull; all 
other bitangent line segments are said to be interior bitangent line segments; 
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the number of boundary bitangent line segments is denoted h; free space is 
the complement in the plane of the interiors of the Oj; a pseudotriangulation is 
a maximal (for the inclusion relation) family of pairwise interior non-crossing 
free bitangent line segments. A pseudotriangulation contains the h boundary 
bitangent line segments plus 3n — 3 — h interior bitangent line segments (thus 
3n — 3 altogether) and induces a decomposition of the free part of the convex 
hull of the Oi into 2n — 2 pseudotriangles [3B] . Figure [3] shows a family of 
7 pairwise disjoint convex bodies of the real affine plane, its (6 in number) 
boundary bitangent line segments, and one of its pseudotriangulation. The 




Figure 3. 

set of pseudotriangulations of a family of convex bodies depends only on its 
chirotope; cf. Appendix [Bl Therefore it is sensible to ask if a pseudotrian- 
gulation of a family of convex bodies presented by its chirotope is efficiently 
computable and, more generally, it is sensible to ask if a pseudotriangulation 
that contains a given set of pairwise interior non-crossing distinguished free 
bitangent line segments is efficiently computable. The main result of the pa- 
per is a positive answer to the first question and, at the same price, a positive 
answer to a restricted version of the second question. 

Theorem 1.1. A pseudotriangulation (and in particular the boundary bitan- 
gent line segments) of a finite planar family of n pairwise disjoint convex 
bodies presented by its chirotope is computable in 0(n log n) time and linear 
space. A similar result holds for the problem of computing a pseudotriangula- 
tion that contains a given set of pairwise interior non-crossing distinguished 
free bitangent line segments, under the assumption that the number of distin- 
guished bitangent line segments that appear consecutively on the boundary of 
any pseudotriangle of any pseudotriangulation of the family of convex bodies 
containing the distinguished bitangent line segments is a constant. □ 

1.1.3. Three independent algorithms. Subsequently we use the term family of 
pairwise disjoint convex bodies with constraints for a finite planar family of 
pairwise disjoint convex bodies together with a, possibly empty, set of pair- 
wise interior non-crossing distinguished free bitangent line segments, the con- 
straints for short; in this context, free space is the space obtained by cutting 
the complement in the plane of the interiors of the convex bodies along the 
constraints: this is the disjoint union of two-dimensional surfaces whose cuffs 
contain exactly one cusp point per endpoint of constraint (counting multi- 
plicities); in particular if the set of constraints is a pseudotriangulation, free 
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space is the disjoint union of the pseudotriangles of the pseudotriangulation 
plus the complement in the plane of the interior of the convex hull of the 
bodies. The family will be said well- constrained if it satisfies the condition 
stated in the theorem above, that is, if the number of constraints that appear 
consecutively on the boundary of any pseudotriangle of any completion of the 
set of constraints into a pseudotriangulation of the family of convex bodies 
is a constant. Our pseudotriangulation algorithm is the composition of three 
independent algorithms: 

(1) an algorithm to compute the convex hull, i.e., the boundary bitangent 
line segments, of a planar family of pairwise disjoint convex bodies; 

(2) an algorithm to compute a cross-section of the visibility complex of a 
family of pairwise disjoint convex bodies with constraints; and 

(3) an algorithm to compute the greedy pseudotriangulation associated to a 
given cross-section of the visibility complex of a family of pairwise dis- 
joint convex bodies with constraints whose set of constraints contains 
the boundary bitangent line segments of the family of bodies. 

Before recalling the definitions of the terms visibility complex, cross-section, 
and greedy pseudotriangulation, we add to our two (non-restrictive) assump- 
tions concerning the boundaries of the convex bodies — recall that one of these 
two assumptions says that the boundaries are free of line segments and the 
other one says that there is exactly one tangent through each boundary point — 
the assumption that the family of convex bodies has no triple tangent. This 
additional assumption is not a restriction on the possible inputs of our al- 
gorithm since for any non-simple chirotope there exists a simple chirotope, 
computable in constant time, such that the non-simple chirotope and the sim- 
ple chirotope have the same set of free bitangent line segments and the same 
set of pseudotriangulations; cf. Appendix [Bl 

1.1.4. Visibility complexes. Let X be a connected component or a union of 
connected components of the free space of a given family of pairwise disjoint 
convex bodies with constraints living in a topological plane A. We denote 
by -C(A) and C^^^A) the spaces of lines and directed lines of A and we take 
for granted that the canonical projection C°^(A) — >• C{A) is a two-covering. 
The space X inherits from the topological point-line incidence geometry of A 
a natural partial topological point-line incidence geometry whose system of 
lines C{X) is defined as the space of pairs {x,£) where i ranges over the space 
of lines of A and x the set of connected components of the pre-image of the 
line i under the canonical projection X — > A, and whose set of incidences 
is the set of point-line pairs (j),{x,l)) € X x >C(X) with p £ x. Note that 
the second component £ of a pair {x,i) € >C(X) is determined by its first 
component x unless x is reduced to a point, which happens precisely when 
X is a cusp point of the boundary of X. Except in the case where X is the 
complement of the interior of the convex hull of the bodies, in which case 
>C(X) is a torus to which is attached, along one of its non trivial closed simple 
curve, a one-punctured disk, the space of lines of X has a natural structure of 
(possibly one-punctured) two-dimensional cell complex: its 1-skeleton is the 
set of tangents to the boundary of X — which includes the lines through the 
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cusp points of the boundary of X — and its 0-skeleton is the set of bitangents 
of X. The visibility complex of X is its space of hnes >C(X) endowed with 
its natural structure of cell complex; furthermore we add to the definition 
that the one-skeleton of the visibility complex is endowed with the orientation 
inherited by duality from the orientation of the underlying topological plane; 
cf . [SSI ESI E] • Similarly we introduce the space (X) of directed lines of X, 
endowed with its natural structure of cell complex together with the natural 
orientation of its one-skeleton inherited by duality from the orientation of 
the underlying topological plane, and we take for granted that the natural 
projection >C°''(X) — ?> £(X) is a two-covering in picture of the two-covering 
C°'^{A) — 7> C{A) and that the cell structure on >C°''(X) is regular contrary, in 
general, to that of C(K). 

Example 1.1. The visibility complex jC(X) of the free space X of a family of 
two disjoint convex bodies Oi,Oj is composed of 

(1) four 0-cells: the four bitangents ti,t2,ts,t4 of the family of bodies; 

(2) eight oriented 1-cells: the four connected components of o*\{ti, ^2, ^3) ^4} 
and the four connected components of o* \ {ti, t2> ^a, ^4}, where o* de- 
notes the set of tangents to Oj ; and 

(3) five 2-cells : the sets of lines with labels — in the context of a family of 
convex bodies with empty set of constraints, the label of a directed line 
is the sequence of bodies intersected by the line ordered as they appear 
along the line and prefixed or postfixed or both prefixed and postfixed 
by the symbol 00 in case the line is (orientation preserving) homeo- 
morphic to M"*", or M endowed with their natural orientations — with 
labels ij, ioo, joo, the set of lines with label cxdoo that separate the 
two bodies, and the set of lines with label 0000 that do not separate 
the two bodies; 

put together as indicated in Figure H] where we write i for the bitangent tf, this 
complex is not regular : the boundaries of the 2-cells with label ioo and joo 
are complete graphs on four elements; and this complex has one end, indicated 
by a marked point 00 (red in pdf color) in the figure. 

Example 1.2. The visibility complex >C(X) of a generic pseudotriangle X with 
cusp points a, b, c consists of 

(1) three 0-cells : the tangents ta,th and tc at the cusp points a,b, and c; 

(2) six oriented 1-cells : the x* = C{x) \ {tx}-, x € {a, 6, c}, where C{x) 
denotes the set of lines through the point x, and the three connected 
components a, /3 and 7 of the curve C{dX) of tangent lines to the 
pseudotriangle minus ta,tb , and tc, and 

(3) three 2-cells : the interiors of the C{x, x') where C{x, x') denotes the 
set of lines joining the sides opposite to the pair of cusp points x and x'; 

put together as indicated in Figure [5l again observe that this complex is non 
regular (its one-skeleton is already non regular); this complex has no end. 

1.1.5. Cross- sections. The boundary of any bounded 2-cell of /^^(X) — bounded 
in the sense that the cell contains no end of >C(X) — has a unique vertex of out- 
degree two and a unique vertex of indegree two; therefore one can speak of 
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the source and sink vertices of a 1- or bounded 2-cell of /^"""(X), and one can 
speak of the left and right boundary chains of a 2-cell. Let C(K)* — > C(K) be 
the inverse image of a universal cover £"(A) of £(A) under the natural pro- 
jection C{X) £(A)— that is C{X)* is the set of pairs {v,l) € C{X) x /:"(A) 
such that the image of v under JC(X) — )■ £(A) coincides with the image of / 
under £"(A) C{A), and £(X)* £(X) is the first projection, cf. [T71 pages 
113-114]— let 0(X) be the set of cells of jCiX)* endowed with the partial order 
generated by the relations 

(1) sour(o") ^ a ^ sink(cr) 

where a ranges over the set of 1- and bounded 2-cells of 0{X) and where 
sour((T) and sink((T) denote the source and the sink of the cell a, let be the 
generator of the automorphism group of the covering C{X)* — >■ >C(X) defined 
by the condition that a -< z^((t): the shift operator for short, and finally let 
J be a maximal antichain of 0(X). The cross-section, denoted r(J), of the 
visibility complex of X at the maximal antichain J is the directed multigraph 
whose set of arcs is the set of 2-cells of J and whose set of nodes is the set of 
0- and 1-cells of J, the source node of an arc being defined as the unique node 
included in its right boundary (if any) and its sink node being defined as the 
unique node included in its left boundary (if any). 

Example 1.3. Let / be a proper filter of the subposet of 0-cells of 0(X). 
Then the set of 1- and 2-cells of 0(X) whose sinks belong to / but not their 
sources is a maximal antichain; the corresponding cross-section is called the 
canonical cross-section associated with the filter /. 

Example 1.4. Figure O depicts a family of 7 convex bodies of the real affine 
plane with one constraint (the bodies are numbered from 1 to 7 and the con- 
straint is the undirected version of the right-right bitangent line segment join- 
ing the third body of the family to the fourth body) and (an upward drawing 
of) the canonical cross-section of its visibility complex associated with the fil- 
ter of the subposet of vertices of 0(X) with angle > 0. The family of convex 
bodies is augmented for each 1-cell e of the cross-section with the horizontal 
line t{e) € e. The horizontal lines t{e) induce a trapezoidal decomposition 
of free space whose trapezoids (23 in number) are in one-to-one correspon- 
dence with the arcs of the cross-section : 21 of these 23 trapezoids are labeled 
in the Figure and these labels are reported on the corresponding arcs of the 
cross-section. 

Example 1.5. Figure [7] depicts a family of 7 convex bodies of the real affine 
plane with one constraint (the bodies are numbered from 1 to 7 and the con- 
straint is the undirected version of the right-left bitangent line segment joining 
the second body of the family to the fourth body) and the canonical cross- 
section associated with the filter of 0-cells of 0(X) generated by the lift in 
C{X)* of the principal filter of any left-left lift in £"(A) of a left-left boundary 
bitangent (the one joining the first body to the second body). The family of 
convex bodies is augmented for each 1-cell e of the cross-section with a line 
t{e) G e. The t{e) induce a trapezoidal decomposition of free space whose 
trapezoids (23 in number) are in one-to-one correspondence with the arcs of 
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the cross-section : 9 of these 23 trapezoids are labeled in the Figure and these 
labels are reported on the corresponding arcs of the cross-section. 
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1.1.6. Greedy pseudotriangulations. One of the key results of the theory of 
visibihty complexes is that the set of sink bitangent line segment^ of the cells 
of a cross-section of a visibility complex is a pseudotriangulation; cf. [31 Theo- 
rem 6, Claim 1]. This pseudotriangulation is called greedy because it can also 
be defined as the set of bitangent line segments of the sequence vi,V2, ■ ■ ■ ,Vh 
of vertices of 0(X) defined inductively by Vi is a ^-minimal element in the 
poset of vertices of the filter generated by the cross-section crossing none of 
the elements of the set {vi,V2, ■ ■ ■ , Vi-i}. 

Example 1.6. Figure [8] depicts the greedy pseudotriangulations associated 
with the two cross-sections introduced in Examples 11.41 and 11.51 Some of 
the labels of the arcs of the cross-sections are reported on the corresponding 
bitangent line segments of the associated greedy pseudotriangulations. 

1.1.7. Declination of the main result. Theorem 11.11 can then be declined as 
follows. 

Theorem 1.2. The convex hull of a planar family of n pairwise disjoint con- 
vex bodies presented by its chirotope is computable (under the guise of the 
circular sequence of boundary bitangent line segments of the family of bodies ) 
in O(nlogn) time and linear space. □ 

Theorem 1.3. The canonical cross-section associated with a given boundary 
bitangent line segment ( as defined in Example \1.5\) of the visibility complex of 
a family of n pairwise disjoint convex bodies with constraints presented by its 
chirotope is computable in 0(n log n) time and linear space. □ 

Theorem 1.4. The greedy pseudotriangulation associated with a given cross- 
section of the visibility complex of a family of n pairwise disjoint convex bodies 
with constraints presented by its chirotope is computable in linear time un- 
der the assumptions that the family is well-constrained and that the set of 
constraints contains the boundary bitangent line segments of the family of bod- 
ies. □ 

Of course it is also sensitive to ask if the (cell structure of the) visibility 
complex of (the free space of) a family of convex bodies with constraints pre- 
sented by its chirotope is efficiently computable. Under the assumption that 
the family is well-constrained, a positive answer to that question is given by 
Angelier and Pocchiola 131 Theorem 1] modulo the efficient computation of 
a cross-section and the efficient computation of its associated greedy pseudo- 
triangulation. (The notion of chirotope used in [3] is finer than the notion 
of chirotope that we are using here — however the algorithmic technique de- 
veloped in [31 page 117], called the xi-Walk procedure, can be adapted to 
the present situation; details on this point will be reported in a different pa- 
per.) Therefore combining our Theorems 11.21 [L3l and 11.41 with Theorem 1 of 
Angelier and Pocchiola [3] we get the following theorem. 



Since there are no tritangent the map that assigns to a free bitangent line segment 
its supporting line realizes a one-to-one and onto correspondence between the set of free 
bitangent line segments and the set of vertices of the visibility complex; thus one can speak 
of the sink bitangent line segment of a 0-, 1-, or 2-cell. 



□ 
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Example 11.41 




Figure 8. 



Theorem 1.5. The visibility complex of a planar family of n pairwise dis- 
joint convex bodies presented by its chirotope is computable in 0{k + nlogn) 
time and linear working space where k is the size of the visibility complex. 
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A similar result holds for the visibility complex of a family of pairwise dis- 
joint convex bodies with constraints under the assumption that the family is 
well- constrained. □ 

In particular the well-constrained chapter of the above result can be used 
to show that the visibility graph of a finite planar family of pairwise interior 
non-crossing line segments presented by the chirotope of the endpoints of the 
line segments is efficiently computable; cf. Appendix O 



1.2. Previous work. The convex hull and pseudotriangulation problems have 
been addressed in the past only for families of pairwise disjoint convex bod- 
ies of an affine topological plane — strictly speaking the problems have only 
been studied in the real affine plane, however it is simple exercise to adapt 
the arguments to afhne topological planes — the following solutions have been 
reported: the set of boundary bitangent line segments can be computed as the 
set of breakpoints of the upper envelope of the support functions of the bodies 
using a divide-and-conquer algorithm, cf. HH chap. 6] and [IQ], and a pseudo- 
triangulation can be computed using a straight sweep a la Bentley-Ottmann 
from the positive horizontal direction to the negative horizontal direction of 
a dynamically changing visibility complex, cf. [37]. Both algorithms run in 
0(n log n) time using not only the chirotope of the family of convex bodies 
but also the direction or slope order on the set of bitangents of the family 
augmented with a point outside the convex hull of the bodies, an information 
which is meaningless in a topological plane which is not affine; the situation 
is even worse for the constrained pseudotriangulation problem since the algo- 
rithm uses also the chirotope of the family of bodies and constraints, that is, 
also the relative positions of the endpoints of the constraints with respect to 
the bitangents. (To fix the ideas we mention that given four pairwise disjoint 
ellipses in the real afhne plane evaluating the position of an endpoint of a bi- 
tangent line segment joining the first two ellipses with respect to a bitangent 
joining the last two ellipses is out of the reach of the current practical tech- 
niques in formal calculus: Grobner bases and so one 06].) More sophisticated 
techniques — using even more involved predicates like slicing the bodies — have 
been developed to design output sensitive convex hull algorithm, cf. [31] • The 
related but different problem of computing the convex hull of a simple curved 
polygon is addressed in [S]. 

We mention that our pseudotriangulation algorithm accepts a larger set 
of inputs, uses simpler data-structures and simpler geometric predicates, has 
fewer degenerate cases to handle, and is faster by a logn factor in its main 
phase (which consists of deriving a pseudotriangulation from a cross-section 
of the visibility complex of the family of convex bodies with constraints) than 
the one developed in |37j and currently implemented in the visibility complex 
package of the CGAL library |2]. 

For families of points the situation is different: Graham's scan [T^ and the 
Knuth's two incremental algorithms [26[ pages 45-61] compute in 0{nlogn) 
time the convex hull of a family of points using only its chirotope; on the 
other hand neither the Chan's output sensitive convex hull algorithm [10] nor 
the one of Kirkpatrick and Seidel |24j are only based on the chirotope since a 
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preliminary step of both algorithms is to compute in linear time an extreme 
point of the family (the one with minimum horizontal coordinate) , a problem 
known to be open for families of points only given by their chirotopes |26[ page 
98] . Similarly a greedy pseudotriangulation of a finite planar family of points 
can be computed in 0(n log n) time using only the chirotope of the family of 
points as we explain in Appendix iDl 

1.3. Outline of our pseudotriangulation algorithm. The design and cor- 
rection of our pseudotriangulation algorithm relies on an extension of the the- 
ory of visibility complexes of families of pairwise disjoint convex bodies of 
the real affine plane to families of pairwise disjoint convex bodies of topolog- 
ical planes and of their branched coverings. In particular our Theorem 11.41 is 
not only valid for families of pairwise disjoint convex bodies with constraints 
of topological planes but also for families of pairwise disjoint convex bodies 
with constraints of branched covering of topological planes (under the mild 
assumption that the convex bodies cover the branch points of the covering 
space). A similar observation can be made regarding Theorem 1 of Angelier 
and Pocchiola [3]. While the use of universal coverings, or portions of universal 
coverings, in the design of geometric algorithms had already appeared in the 
early days of the computational geometry literature, e.g., it seems to 

be the first time that branched coverings are used in the design of a geometric 
algorithm. (Branched coverings are used in |43] to define the dual Voronoi 
diagram of a constrained Delaunay triangulation in the plane, but apparently 
without algorithmic consequences — see also the discussion in [T21 page 30].) 
We refer to [28l page 145], [23], [271 page 18] and the references cited therein 
for background material on branched coverings. 

Our algorithm proceeds in three steps: we first compute the convex hull of 
the family of convex bodies, then the cross-section of the visibility complex 
of the family of convex bodies with constraints assigned to a distinguished 
boundary bitangent line segment, and finally the greedy pseudotriangulation 
associated with that cross-section, that is, the set of sinks of its 2-cells, cf. \37\ 
Theorem 12] and more generally [31 Theorem 6, Claim 1] in the case where we 
look for a constrained pseudotriangulation. 

1.3.1. Convex hull algorithm. Our convex-hull algorithm is a sweep of a con- 
nected 4-sheeted branched covering of the underlying plane ramified over any 
interior point of an arbitrarily distinguished body: we sweep the 4-sheeted 
covering surface with a half-line whose supporting line is a left tangent at the 
origin of the half-line to the lift of the distinguished body. Any body, except 
the distinguished one, has four lifts in the 4-sheeted covering surface; we only 
keep the lifts either lying in one of the first three sheets, either straddling 
the first two sheets or the second and third sheets or the last two sheets, as 
illustrated in Figure [9] where the bodies numbered 1, 3, 6 and 7 are lifted only 
in the first three sheets and where the bodies numbered 2 and 5 are lifted 
astride the first two sheets, the second and third sheets, the last two sheets 
but not astride the last and first sheets. The sweep starts at a boundary tan- 
gent and induces a total order on the lifted bodies with the property that a 
body contributes to one or zero connected piece to the boundary of its convex 
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Figure 9. 

hull with its predecessors in the total order. During the sweep we maintain the 
convex hull of the lifts of the bodies that have been entirely swept or partially 
swept by the sweeping half-line; the convex hull of the family of bodies is then 
extracted from the convex hull of the lifts, as illustrated in Figure [9] where one 
can read the convex hull of the family of bodies as the boundary bitangents 
of the lifts drawn with a bold line. 

1.3.2. Cross-section algorithm. Our cross-section algorithm is again a sweep 
but now a simple sweep of the convex hull of the bodies by a half-line whose 
supporting line is a left tangent at the origin of the half-line to one of the bodies 
appearing on the boundary of the convex hull — a boundary body, for short. 
The sweep starts at one of the boundary bitangent line segments leaving the 
distinguished boundary body. During the sweep we construct the canonical 
cross-section of the visibility complex of the family of convex bodies with 
constraints assigned to the distinguished boundary bitangent line segment, 
cf. Example II. 5[ The method presents some interesting and novel features 
due to the fact that the relative positions of the constraints with respect to 
the bitangents are not completely determined by the chirotope of the convex 
bodies. It is also interesting to mention that this second step is implementable 
in O(nlogn) without restriction on the possible sets of constraints. 

1.3.3. Greedy pseudotriangulation algorithm. Our third and last algorithm — 
which consists of deriving the greedy pseudotriangulation associated to a given 
cross-section of the visibility complex of a family of pairwise disjoint convex 
bodies with constraints whose set of constraints contains the boundary bi- 
tangent line segments of the bodies — is the most elaborate and fully benefits 
from the idea of using branched coverings. A preliminary version of this third 



16 



LUC HABERT AND MICHEL POCCHIOLA 



algorithm — of which the idea of using branched coverings was unfortunately 
missing — was discussed several years ago by the second author of the pa- 
per with his PhD student Pierre Angelier, see [H pages 83-92] and compare 
with [33 Appendix A]. 

We define a partial order < on the set of 2-cells a of the input cross-section 
whose sink bitangent line segment t{a) is not a constraint (and thus not a 
boundary bitangent line segment), and for each a we define a pair of adjacent 
pseudotriangles, called the ^;S-pseudotriangles of a, made with the t(cr'), a' < 
a, and with auxiliary bitangent line segments Sj{a), 1 < J < cr*, such that 
a representation of the ^,B-pseudotriangles of cr by a linked structure TZa — 
that is, collections of nodes interconnected by pointers; cf. |45l page 8] — is 
computable in constant amortized time and such that t{a) is computable as 
the bitangent line segment joining the ^;B-pseudotriangles of a in constant 
amortized time starting from the knowledge of the linked structure TZ^- 

A key feature of our method is that the ^S-pseudotriangles are defined as 
projections in the plane of pseudotriangles of pseudotriangulations of sets of 
lifts of bodies in certain branched coverings of the plane. (Some of the Sj{a), 
1 < i < are computed by a recursive application of the procedure to 
compute the t{a).) More precisely, given a finite family of pairwise disjoint 
convex bodies with constraints (including the boundary bitangent line seg- 
ments of the convex bodies) of a branched covering B of a topological plane 
A, we associate to each bounded 2-cell a of its visibility complex whose source 
bitangent line segment is not a constraint a pseudoquadrangle containing |J a, 
called the 'H-pseudoquadrangle of a and denoted T-L^a), whose diagonals are 
the source and the sink bitangent line segments of a; pseudoquadrangle from 
which we derive, once a cross-section T containing a is chosen, a pair of pseu- 
dotriangles adjacent along the source bitangent line segment of a, called the 
^i3o-pseudotriangles of a, with the property that the bitangent line segment 
joining the ^^So-pseudotriangles of a is the sink bitangent line segment of a; 
the definition of the >i,Bo-pseudotriangles depends on the type of <t in P which 
is a pair ij, i,j G {1,2,3}, that encodes the position of the source and sink 
nodes of a in the decomposition of the left and right boundaries of a into 
convex chains (3 in number at most). Then we assign to the 2-cell cr, element 
of the cross-section P, a 2-cell /^(ct), element of a certain cross-section /u(P) of 
the visibility complex of a certain family of convex bodies and constraints of a 
certain branched covering /i(B) of the topological plane A — obtained as con- 
nected sum of B and copies of the plane A as indicated in Figure [10] — so that, 
among other things, a and /^(ct) have the same sink. The ^;S-pseudotriangles 
of a are then defined as the ^0o-pseudotriangles of n{a). The correction 
of the method relies on several new properties of cross-sections of visibility 
complexes. 

1.4. Organization of the paper. In the next section we extend the theory 
of pseudotriangulations and visibility complexes to the setting of branched 
coverings of topological planes (no proofs will be given since one can adapt 
easily to that setting the proofs given in [381 133 [3] ) , we establish several new 
properties of cross-sections of visibility complexes, and we introduce the main 
ingredients of our pseudotriangulation algorithm mentioned in the previous 
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Figure 10. Given a branched covering B of a topological plane A 
and a simple oriented curve 7 in B homeomorphic to its projection 
7' in A under the covering map B A we construct a new branched 
covering A^(B) of A as follows. Cut C = BUA along 7 and 7', call C-,, 
the resulting surface and q : <C-y — > C the induced projection. Then 
we define A^(B) as the quotient space of by identification of the 
left lift of 7 under q with the right lift of 7' under q as well as the 
right lift of 7 under q with the left lift of 7' under q. The operator /i 
is defined as the composition of several operators. 

sections. In the third section we describe our pseudotriangulation algorithm, 
we analyze its complexity, and we conclude in the fourth and last section. 



18 



LUC HABERT AND MICHEL POCCHIOLA 



2. Visibility in branched coverings 

In this section we extend the theory of pseudotriangulations and visibihty 
complexes to the setting of branched coverings of topological planes; we es- 
tablish several new properties of cross-sections of visibility complexes; and we 
introduce the key ingredients of our algorithm mentioned in the introduction : 
?^-pseudoquadr angles, ^i3o-pseudotriangles, and ^S-pseudotriangles. For the 
sake of simplicity and clarity we only went over the case where the set of 
constraints is empty, the general case can be treated very similarly using the 
definition of visibility complexes of families of pairwise disjoint convex bodies 
with constraints given in the introduction. 

Let V he a finite family of pairwise disjoint convex bodies of a finite con- 
nected branched covering space B of an oriented topological plane A equipped 
with the partial topological point-line incidence structure, with singularities at 
the branch points, inherited from the point-line incidence structure of A. We 
assume that the boundaries of the convex bodies are free of line segments, that 
there is exactly one tangent line through each boundary point, that the bodies 
surround the branch points of the covering space, and we use the following 
associated terminology and notations: free space is the complement of the in- 
teriors of the bodies; a bitangent line segment is a closed line segment of free 
space tangent to two bodies at its endpoints; a boundary bitangent line segment 
is a bitangent line segment contained in the boundary of the convex hull of 
the bodies; all other bitangent line segments are said to be interior bitangent 
line segments; a primitive arc is a connected component of the boundary of 
the bodies minus the bitangent line segments; hx) is the number of boundary 
bitangent line segments; nj) is the sum of the orders of the branch points plus 
the number of bodies surrounding no branch points; k-z) is the number of sheets 
of the branched covering space. 

2.1. Pseudotriangulations. A pseudotriangulation is a maximal, for the in- 
clusion relation, collection of pairwise interior non-crossing bitangent line seg- 
ments. As in the case where the covering map IB ^ A is the identity map of 
the real affine plane, a pseudotriangulation induces a subdivision of free space 
whose bounded regions are pseudotriangles, that is, subsets of free space home- 
omorphic via the covering map to pseudotriangles of the topological plane. 

Theorem 2.1. Let T be a pseudotriangulation ofD. Then the bounded faces 
of the subdivision of free space induced by T are pseudotriangles, their number 
is 2nx) — 2kx> and the size of T is Sn© — Sfcxi- Furthermore any interior 
bitangent line segment of T can be flipped, that is, replaced by an interior 
bitangent line segment to obtain a new pseudotriangulation. □ 

Proof One can repeat the proof given for the real affine plane in [38j since 
the lines of a topological plane — and consequently the lines of free space — are 
geodesies for an ad hoc metric on the topological plane; cf. [9l Theorem 11.2, 
page 56]. □ 

Two pseudotriangulations are said to be adjacent (or related by a flip) 
if they differ by a single (necessarily interior) bitangent line segment. The 
adjacency graph on the set of pseudotriangulations is a connected regular 
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graph of degree 3nx> — Skv — hx>. More generally the collection, ordered by 
inclusion, of subsets of pairwise interior non-crossing free interior bitangent 
line segments is a strongly flag-connected pure simplicial complex of dimension 
2>nx> — 3kxi — hx> which satisfies the diamond property. This simplicial complex 
will be called thereafter the complex of pseudotriangulations of the family of 
convex bodies. 

Example 2.1. Figure [TT] depicts a family of two convex bodies of a 2-sheeted 
branched covering of A with two branch points (the two sheets are obtained by 
cutting the covering space along the two line segments joining the two branch 
points). Its complex of pseudotriangulations is the cocube of dimension 2. 




Figure 11. 

More generally the complex of pseudotriangulations of a family of two convex 
bodies of a n-sheeted covering surface of the plane with two branch points is 
the cocube of dimension n. 

2.2. Visibility complexes. We now assume that there is no tritangent. Free 
space in denoted F. The space F inherits from the point-line incidence struc- 
ture of A a natural partial point-line incidence structure whose system of lines 
is defined as the space of connected components of the pre-images of the lines 
of A under the canonical projection X ^ A, and whose set of incidences is the 
set of point-line pairs (p, x) with p G x. The label of a directed line of F is the 
sequence of bodies intersected by the line ordered as they appear along the 
line and prefixed or postfixed or both prefixed and postfixed with the symbol 
oo in case the line is (orientation preserving) homeomorphic to the curves M^, 
M~ or R endowed with their natural orientations. A directed line of F touch- 
ing tangentially a body o is called a left or right tangent to o depending on 
whether o lies, locally around the touching point, on the left side or on the 
right side of the line. A directed line of F joining tangentially a body o to a 
body o' is said to leave o and to reach (or enter) o' and is called a left-left, 
left-right, right-left, or right-right bitangent depending on whether the line is 
a left tangent to both o and o', a left tangent to o and a right tangent to o', 
a right tangent to o and a left tangent to o', or a right tangent to both o and 
o' . The sets of left and right tangents to a body are simple closed curves to 
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which we assign the orientation inherited by duaUty from the orientation of 
the ground topological plane A. 

2.2.1. Cell structure. Let V = be the space of directed lines of F, V"^ its 
space of left and right tangents, and V'^ its space of left-left, left-right, right- 
left and right-right bitangents. The operator that reverses the direction of a 
directed line is denoted l and we take for granted that the natural projection 
V ^> V/t is a 2-covering. The increasing sequence 

(2) c C = V 

is, modulo the adjunction of a point at infinity in each connected component 
of V\V^ whose topological closure is noncompact, the sequence of 0-, 1-, and 
2-skeletons of a natural structure of finite 2-dimensional regular cell complex 
on Y: since the curves of tangents to the bodies are oriented curves one can 
speak of the source and sink vertices or 0-cells of a 1-cell; as usual a chain 
of is a sequence of 0- and 1-cells such that the predecessor (if any) and 
the successor (if any) of a 1-cell are its source and its sink, respectively; as 
usual the points added at infinity are called the ends of V; and a 2-cell is said 
bounded if it contains no end. This complex satisfies the following properties: 

(1) A 0-cell is the source and the sink of two 1-cells; 

(2) The boundary of a bounded 2-cell is composed of two chains that share 
the same source/sink, called the source/sink of the 2-cell. Conversely 
any vertex is the source/sink of a bounded 2-cell. By convention the 
right /left boundary chain of a bounded 2-cell a with source v, denoted 
rc(<T) /1c((t), is the boundary chain of a whose first 1-cell is supported by 
the curve of tangents to the body reached/left by v and supporting v; 

(3) The boundary of an unbounded 2-cell is composed of a single chain; 
An unbounded 2-cell is said to be left or right unbounded depending 
on whether its boundary is composed of right or left tangents, respec- 
tively. The number of left unbounded 2-cells and the number of right 
unbounded 2-cells are both equal to the number of connected compo- 
nents of the complement of the convex hull of the family of convex 
bodies; 

(4) a 1-cell is incident to three 2-cells with labels the subsequences of length 
two of its label; The three 2-cells incident to the 1-cell a of V with label 
ijk are denoted (Tq, a G {r, 1, b, f}, according to the following rule: the 
2-cells with label ij and jk are denoted o"b and o"f, respectively; and 
the remaining 2-cell, whose label is ik, is denoted cJr or ai depending 
on whether the lines of a are left or right tangents, as illustrated in 
Figure ^ 

(5) a 0-cell is incident to four 1-cells with labels the subsequences of 
length 3 of its label. 

(6) a 0-cell is incident to six 2-cells with labels the subsequences of length 2 
of its label. The six 2-cells incident to the 0-cell a of V with label ijkl 
are denoted (Jq, a G A = {e, s, r, 1, b, f} according to the following rules: 

(a) (Tb is the 2-cell with label ij; 

(b) (Tf is the 2-cell with label kl; 
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Figure 12. 



(c) is the 2-cell with label jk,jl,il, or ik depending on whether a 
is a left-right, left-left, right-left or right-right bitangent; 

(d) (Tg is the 2-cell with label il,ik,jk, or jl depending on whether a 
is a left-right, left-left, right-left or right-right bitangent; 

(e) (Tr is the 2-cell with label ik,il,jl, or jk depending on whether a 
is a left-right, left-left, right-left or right-right bitangent; 

(f) 0"! is the 2-cell with label jl,jk,ik, or il depending on whether 
£7 is a left-right, left-left, right-left or right-right bitangent, as 
illustrated in Figure [13l 




Figure 13. 



By definition the visibility complex of the family of pairwise disjoint convex 
bodies T> is the regular cell-complex V'^ c V"^ C = V endowed with the 
orientation of its one-skeleton inherited from the orientation of the ground 
topological plane A. 

Example 2.2. The visibility complex a family of two disjoint convex bodies 
Oj, Oj of A is composed of eight 0-cells, sixteen (oriented) 1-cells and ten 2-cells 
(the sets of lines with labels ij, ji, ioo, ooi, joo, ooj, and oooo four times) put 
together as indicated in Figure [TH where, by convention, the left boundary 
chain of a bounded 2-cell is above its right boundary chain (thus, one can 
read on the Figure that the right boundary chain of the 2-cell with label ij 
is lji3'i'^4:' and that its left boundary chain is Iii2j24'); in the introduction 
section we observed that the quotient of this complex under t is not regular. 
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Figure 14. 

Example 2.3. Figure [TSl depicts the cell decomposition of the quotient under 
t of the visibihty complex of a family of two convex bodies of a 2-sheeted 
branched covering of A with two branch points (the two sheets are obtained 
by cutting the covering space along the two hne segments joining the two 
branch points). 



January 9, 2013 



1 ii-^ 3 2 




1 

1 


Jl ^ 


7 


7 




— 6 


6 




^ 4 


4 




^ 5 


5 


c- is^— 


^ 3 


3 




^5 2 


2 




^ 8 


8 




^ 1 




Figure 15. 



24 



LUC HABERT AND MICHEL POCCHIOLA 



2.2.2. Horizon operators. We now describe, in preparation for the section on 
^;B-pseudotriangles, the boundary chains of the 1- and 2-cens in terms of the 
operators, denoted ipa, a £ {e, s, r, 1, b, f}, that assign to a 0- or l-ceh a the 
sink vertices of its incident 2-cens aa (if defined); in particular, ipsiv) is the 
sink of the 2-cen with source v, (peiv) is the identity operator, and, for a left 
1-cell e, supported by the curve of left tangents to the body o, the bitangent 
(ff (e) is the first bitangent leaving o encountered when we traverse the curve of 
left tangents to o starting from e. Note that fi and are the conjugates of (/?r 
and under the reorientation operator l, that is, LOip^ = ip^-oL, ioip^ = (p]~,oi. 
We name these operators the horizon operators in reference to the operators 
underlying the definition of the horizon trees of Edelsbrunner and Guibas |T3] . 
For example the table of the horizon operators on the set of bitangents of the 
visibility complex of two convex bodies of the plane is the following 



where - stands for undefined and where we use the notations of Example 12.21 
The proofs of the two following theorems are easy (using continuity arguments) 
and are left to the reader. 

Theorem 2.2. Let e be a left 1-cell supported by the curve of left tangents 
to the body o. Then sink(e) = y'r(e) i/sink(e) reaches o; otherwise sink(e) = 
(/jf (e). Furthermore (p{{e) is the first 0-cell leaving a encountered when we tra- 
verse its curve of left tangents starting from e. A similar result holds for right 
1-cells using conjugation under t. □ 

Theorem 2.3. Let a be a bounded 2-cell of V, let a be the body that the 
source of a reaches, let a' be the body that the sink of a leaves, let c be the 
curve of tangents to a supporting the source of a, and let d be the curve of 
tangents to a' supporting the sink of a. Then the right boundary chain of a 
is the concatenation of three (convex) chains rci(cj), rc2(cr) and rc3{a) whose 
atoms a, except sour((T) and smk{a), are characterized by (ph{a), (pi{a) and 
iff {a) = sink(cr); respectively. Furthermore 

(1) if c is the curve of right tangents to the body o then rci{a) = sour((T); 
otherwise 

rci((T) = sour(o-)eio-uiieii . . . eik^, (fei > 0), 

where ^11^12 • • • viki is the maximal sequence of consecutive 0-cells leav- 
ing o that follow sour(cr) on c; 

(2) if c' is the curve of right tangents to the body d then rcs (cr) = sink((j); 
otherwise 



1 2 3 4 1' 2' 3' 4' 



1' 2' 3' 4' 1 2 3 4 

4' 3' 2' 1 4 3 2 1' 

2 - - 2' 3 4 4' 3' 
3' 4 1' 3 2' - - 2 
2' 2' 2 2 3' 3' 2' 2' 

3 3 2 2 2 2 2' 2' 



rc3(cj) = e3fc3 . . . esit-sieso sink(o-), (A;3 > 0), 
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where v^ks ■ ■ ■ V31 is the maximal sequence of consecutive 0-cells reach- 
ing d that "precede sink((j) on d . 
(3) if c = c' then c is the curve of right tangents to o and 102(0") is a 1-cell 
whose source and sink are the source and the sink of a, respectively; 
otherwise 

rc2{a) = 620^21621^22 • ■■V2k2e2k2, {h > 1), 

where 

(a) 620 is the empty chain if and only if c is a curve of left tangents; 

(b) 62^2 is the empty chain if and only if c' is a curve of left tangents; 

(c) V21 is the first 0-cell reaching that follows sox\x{(t) on c; 

(d) V2k2 is the first 0-cell leaving o' that precedes sink(cr) on c' ; 

(e) ^2,1+1 = Vb(f^2i) o.nd the e2i are right 1-cells. 

A similar result holds for the left boundary chain of a using conjugation un- 
der L. □ 



Example 2.4. The convex decompositions of the left and right boundary 
chains of the bounded 2-cells of the visibility complex of two convex bodies of 
the plane are given in the following table 



a 


rci((j) 


rC2 (cr) 


rc3(o-) 


1ci((t) 


lC2(a) 


lC3(^) 


ij 


Iji 


3' 


^',A' 


lii 


2 


i24' 


ji 


l'i[ 


2' 


J24 




3 


«34 


ioo 


2 


J24' 


^4 ^ 2 ^2 ^ 


2i23i34i4 


Iji 


3' 


ooi 


2'i'23'i'^A'i'^ 




3 


2' 


J^4 


i4lii2i23 


joo 


3 




j4lji3'i^2' 


3j32j24'j^ 


I'ii 


2' 


ooj 


3'j^2'j^4j4 


K 


2 


3' 




jU'jl3j32 


0000 


4' 




1' 


4' 




1' 


0000 


4 




1 


4 


k 


1 



where we use the notations of Example 12.21 



Example 2.5. Consider the family of 7 convex bodies 01,02, ... ,07 of the 
real affine plane depicted in Figure [16] and let a be the 2-cell of its visibility 
complex that contains the directed line labeled cr. Then, using the notations 
bij,bi-,b--7,b77 for the left-left, right-left, left-right, right-right bitangents join- 
ing Oi to Oj, its source and sink are the bitangents 621 and 6yj and the convex 




Figure 16. 
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decompositions of its left and right boundary chains are given in the foUowing 
table 



rci(fT) rc2((T) 


rc3(o-) 


lci(a) 


1c2(ct) 


1c3(ct) 


^2lftl7^17&16&i6^15^14 %i&76 


hi 


&21 


623% 


%4 ^3 ^^2 



where we only indicate the bitangents of the chains. 



2.2.3. Greedy pseudotriangulations. Let L be the set of directed lines of A. 
Let W — ^ V be the inverse image of a universal cover L of L under the natural 
projection V — )■ L — that is, W is the set of pairs {v,l) G V x L such that the 
image of v under V ^ L coincides with the image of / under L ^ L, and 
W — 7> V is the first projection, cf |T71 pages 113-114] — let O be the set of cells 
of W endowed with the partial order generated by the relations 

(3) sour(o") ^ a ^ sink(cr) 

where a ranges over the set of 1- and bounded 2-cells of O and where sour((T) 
and sink((T) stand, respectively, for the source and the sink of a. The sets 
of left and right unbounded 2-cells of W are denoted and 1, respectively; 
note that the elements of and 1 are isolated elements in O and that the 
sizes of and 1 are both equal to twice the number of sheets of the branched 
covering space B. Finally we denote by i' the generator of the (infinite cyclic) 
automorphism group of the covering W — > V/i defined by the condition that 
a -< the shift operator for short, and we keep the same symbol to denote 
a horizon operator and its lift in W; thus ips is the map that assigns to a vertex 
V of O the sink of the 2-cell of O whose source is v. Two vertices of O are said 
crossing if their corresponding bitangent line segments are crossing. 

Theorem 2.4 ([3l Theorem 5] and [37\ Lemma 8]). Two crossing vertices are 
comparable with respect to the partial order ~< and the map 93 : — )■ that 
associates with v G the minimum element of the set of u G such that u 
crosses v and v < u is well-defined, one-to-one and onto. Furthermore if v is 
an interior vertex then ip{v) = (ps{v); otherwise <f{v) = ^{v). □ 

Let J be a maximal antichain of O, let J'^ be the filter of cells 2 € O such 
that X ^ z for some x G J, and let 

(4) G{J) = [jB,{J) 

i>l 

where Bi{J) is the set of minimal elements of the set of vertices of J"*" that do 
not cross any element of |Ji~^ Bj{J) where as usual Ij'j' Bj{J) = 0. We denote 
by b the operator that assigns to a vertex of O its corresponding bitangent 
line segment. 

Theorem 2.5 ([3l Theorem 5] and [371 Theorem 12]). Let J be a maximal 
antichain of O. ThenhoG{J) is a well-defined pseudotriangulation and 

(5) G(J) = 4 \ ip{J+) = V( J) + sink(7-°(J)) 

where Jq is the set of vertices of , V(J) is the set ofO-cells of J, andT°{J) 
is the set of bounded 2-cells of J minus the Cf and cr'^ where cr ranges over the 
set of right-right boundary 0-cells and right boundary 1-cells of J and where 
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a' ranges over the set of left-left boundary O-cells and left boundary 1-cells 
of J. □ 

The pseudotriangulation G( J) is called the greedy pseudotriangulation at J. 

We describe, again in preparation for the section on ^;B-pseudotriangles, 
the boundary chains of the pseudotriangles of the greedy pseudotriangulations 
in terms of the horizon operators. Let J be a maximal antichain of O. Let v 
be a minimal element of the subposet of vertices of J"*" that is not a left-left 
boundary bitangent and let R(f ) be the pseudotriangle of the pseudotriangula- 
tion b o G( J) lying locally on the right side of the bitangent line segment h{v). 
One can easily show that the pseudotriangle R(f ) is independent of the choice 
of the maximal antichain J. Walking in counterclockwise order along the 
boundary of R(f ) starting at the tail of h{v) we traverse successively 4 convex 
chains R-' (v) [j = 1, 2, 3, 4). A description of these chains in terms of horizon 
operators is given in the following theorem where Tp^ denotes the conjugate of 
the operator under b . 

Theorem 2.6 ([3l Theorem 10]). Let v be a bitangent that is not a left-left 
boundary bitangent and let 

BJ{v) = ejoVjieji, . . . , Vju^ejUp kj > 

where Vjk stands for a bitangent line segment and ejk for an arc. Then 

(1) vii = h{v) and vij+i = ^f(fij); 

(2) f3i = Tpy^{vii) and vs^i+i = Tpf{v3i) (assuming that V31 is well-defined); 

(3) V21 =^i.(uii) and V2,i+i =^b(^2i) (assuming that V21 is well-defined); 

(4) Tp{v) leaves an arc o/R^(v) or the first arc ofR^{v). □ 

2.2.4. Cross-sections. Let J be a maximal antichain of O and let V(J), £{J), 
J^(J), and Vl{J) be its sets of 0-, 1-, bounded and unbounded 2-cells. Using 
the simple fact that a maximal antichain and a maximal chain intersect in a 
single element one can easily check that 

(1) U{J) is the whole set of unbounded 2-cells and its size is 2kx>. 

(2) for any bounded 2-cell o" of J there is exactly one atom (a 0- or 1-cell) of 
its right/left boundary chain — denoted botj(o")/topj(o") thereafter — 
that belongs to J; 

(3) for any left/right unbounded 2-cell cr of J there is exactly one atom (a 
0- or 1-cell) of its right/left boundary chain — denoted bot j{a) /top j{a) 
thereafter — that belongs to J; 

(4) the 2-cells of J are exactly the cir, Cf , (Ti, and Cb, where a ranges over 
V( J) + £"( J) and where by convention we ignore or (exclusive) ai if 
one of them is not defined, that is, if o" is a right 1-cell or a left 1-cell; 

(5) the size of 8{J) is 2nx> - 2#V(J); 

(6) the size of T{J) is 3nx> — kjy — #V{J). 

The cross-section of the visibility complex of the family of convex bodies V 
at the maximal antichain J, denoted r(J), is the directed multigraph whose 
set of nodes is the set of 0- and 1-cells of J and whose set of arcs is the 
set of 2-cells a of J directed from hotj{a) to topj(cr). We use the notation 
Tij{J), i,j G {1,2,3}, for the set of cr G J^iJ) such that botj(cr) S rcj(o"), 
top J (cr) € lcj(cj), and both botj(cr) and top j(cr) are 1-cells; the pair ij is 
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called the type of the arc a; the type of an arc captures exactly the upward 




Figure 17. 

embedding in the plane of this arc together with its adjacent arcs with the 
property that the arcs incident to a node a appear in circular order Ur, cif , <ti, 
and (Tb where by convention we ignore or (exclusive) a\ if one of them 
is not defined, that is, if a is a right 1-cell or a left 1-cell, as illustrated in 
Figure [T71 We make the set of cross-sections into a poset A{0) by defining 
r(J) ^ r( J') in A{0) by J+ ^ J'+ where J+ is the filter of cells z eO such 
that X ^ z for some 2; € J. The covering relations in A{0) are described in 
the following theorem, from which it follows by induction, starting from the 
obviously acyclic cross-sections of Examples 11.41 and 11.51 that cross-sections 
are acyclic. 

Theorem 2.7. Let J be a maximal antichain of O, let v € , e, e', /, f £ , 
a, a' S with v = sink((T) = sink(e) = sink(e') = sour((T') = sour(/) = 
sour(/') . Then v is minimal in the subposet of vertices of J"*" if and only if 
either (first case) a,e,e' £ J or (second case) v £ J. Furthermore in the first 
case J' = J — {fj, e, e'} + v is a maximal antichain and r(J') covers r(J); 
in the second case J' = J — v + {/, f'jCr'} is a maximal antichain and T{J') 
covers T{J); and in both case r(J') is obtained from T{J) by local changes 
as indicated in Figure [721 where the arcs numbered 1, 2, 3, 4 stand for the four 
2-cells Vr,V{,v\, and V]^ incident to vertex v and where the arcs are oriented 
upward. □ 

Theorem 2.8. Let J be a maximal antichain of O. Then T{J) is acyclic 
( with set of sources and set of sinks 1 ) and can be embedded in free space in 
such way that the arcs incident to a node a appear in circular order cjr,cjf,cri, 
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Figure 18. 

and (Tb where by convention we ignore fii- or (exclusive) ai if one of them is 
not defined, that is, if a is a right 1-cell or a left 1-cell. □ 

Example 2.6. The Hasse diagram of the poset of cross-sections of the visi- 
bility complex of a family of two convex bodies is depicted in Figure [191 the 
diagram is of course invariant under the shift operator and its quotient modulo 
the shift operator is composed of twelve cross-sections. 
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Figure 19. The Hasse diagram of the poset of cross-sections of the visibihty complex of a family of two bodies: the arcs are 
oriented upward and the ones incident to a node a are arranged in the circular order w^, Wf , ui, Ub, • • ■; the nodes with degree 3 
are 1-cells; the nodes with degree 4 are 0-cells, that is, bitangents; the arcs incident to a unique node are the unbounded 
2-cells. 
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2.3. ^,B-pseudotriangles. The definition of the ^jB-pseudotriangles is given 
in terms of the horizon operators (/J^ and their dual horizon operators ipa, which 
are defined in exactly the same way except that we replace in the definition of 
the horizon operators the sink operator by the source operator, i.e., (pa,{a) is 
the source of the 2-cell a a', note that ifs^a, = ^a- More precisely we are going 
to use the derived horizon operators Pforw and Pback (and their duals) which 
are defined as follows: 

(1) Pforw is the operator that assigns to a bitangent line segment v its 
image under Tpf if t) is a right-left or left-left bitangent line segment 
whose image under Tp^ is a not a boundary bitangent line segment; the 
bitangent line segment v otherwise; and 

(2) Pback is the operator that assigns to a bitangent line segment v its 
image under Ip-^^ if u is a right-right or left-right bitangent line segment 
whose image under Tp^ is a not a boundary bitangent line segment; 
the bitangent line segment v otherwise; note that pback and piorw are 
conjugate under the reorientation operator l. 

As a simple consequence of the description of the greedy pseudotriangulations 
in terms of the horizon operators given in Theorem 12.61 we see that greedy 
pseudotriangulations are stable under the operators pforw and Pback- 

2.3.1. T-L-pseudoquadrangles. For a a bounded 2-cell we set Hrc(<7) = h{u) if 
rc2(cr) is reduced to a 0-cell u; otherwise we choose a 1-cell e of rc2(<T), we 
introduce the sequence uq = b ip\^{e) , Ui+i = Phackiui), the sequence vq = 
hipf^{e),Vi+i =Pforw.(i'i), and we set 



where H^((T, e) is, depending on whether uq is a boundary bitangent line seg- 
ment or not, the empty sequence or the sequence of Ui truncated just after the 
first index i such that Ui = Uj+i, and where, similarly, H~(o", e) is, depending 
on whether vq is a boundary bitangent line segment or not, the empty sequence 
or the reversal of the sequence of Vi truncated just after the first index i such 
that Vi = Uj+i. The bitangent line segments corresponding to the vertices of 
rc2((T) are consecutive elements of the sequence Hrc(o"), from which it follows 
that Hi.c(<7) is independent of the choice of the 1-cell e. It is convenient to 
extend the definitions of H~(o", .) and H^(fT, .) to the whole set of 1-cells of 
rc((T) as follows 



where e is the empty chain, so that equation ([6]) holds for any 1-cell e of 
rc((7). Let Uq, Ui, . . . be the sequence of bitangent line segments defined by uq 
is the source bitangent line segment of a and tij+i = Pforw«(^i); similarly, let 
vo,vi, . . . be the sequence of bitangent line segments defined by vq is the sink 



(6) 



Hrc((T) = Hrc(fT,e' 



:)H+(a,e) 




Hrc(o") if e is a 1-cell of rc3(fT); 
e if e is a 1-cell of rci (cr); 




Hrc(o") if e is a 1-cell of rci (cr); 
e if e is a 1-cell of rc3(fT); 
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bitangent line segment of a and Vi^i = Pback(^^i); clearly by construction one 
has 

(7) Hrc(a)=Hj,(a)H?,(a)H?,(a) 

where (1) HJj,((t) is the empty sequence if ui is a boundary bitangent line 
segment or if uq = ui; the reversal of the sequence ui,U2,- ■ ■, truncated just 
after the first index i such that Ui = Uj+i, otherwise; (2) H^c('^) is the sequence 
of bitangent line segments corresponding to the sequence of vertices of the 
chain rc2 (cr); and (3) H^p(cj) is the empty sequence if fi is a boundary bitangent 
line segment or if vq = vi; the sequence vi,V2, ■ ■ ■, truncated just after the first 
index such that Vi = fi+i, otherwise. We consider Hj-c as an operator on the 
set of bounded 2-cells and we define Hic to be its conjugate under t. 

Example 2.7. The table of the operators H*^, on the 2-cells of the visibility 
complex of two convex bodies of the plane is the following 
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where we use the notations of Figure [TJ] and where e stands for the empty 
sequence. 

Example 2.8. Consider the family of 7 convex bodies oi, 02, 03, . . . , 07 of the 
real affine plane depicted in the left part of Figure [20] and let a be the 2-cell of 
its visibility complex that contains the directed line segment labeled a. Then, 




Figure 20. 

using the notations ^ii, fc,-, fer for the left-left, right-left, left-right, right- 
right bitangent line segments joining Oj to Oj, its source and sink bitangent 
line segments are the line segments 626 and 643 and the Hl^{a) and HJ^(cr) are 
given in the following table 



Hi,(a) H?,(a) H?,(a) 








^2 ^23 


^34 
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Observe, as illustrated in the right part of the figure, that the bitangent line 
segments of Hrc(o") and Hic(o") are the bitangent line segments of a pseudo- 
quadrangle with diagonals the source and sink bitangent line segments of a. 

Example 2.9. Consider the family of 7 convex bodies 01,02,..., 07 of the 
real affine plane depicted in Figure [21] and let a be the 2-cell of its visibility 
complex that contains the directed line segment labeled a. Then its source 




Figure 21. 

and sink bitangent line segments are the boundary bitangent line segments 621 
and hji and the H*j,(cj) and H{^(o") are given in the following table 



Hj,(a) Ulia) Ulia) 








e 


^23% e- 



Observe that the bitangent line segments of Hrc(o")/ Hic((t) augmented with 
the sink/source of a bound a pseudotriangle. 

We now reinterpret the sequence Hi.c(c) and Hic((t) through the greedy 
pseudotriangulations. Let u be a bounded 2-cell, let e be the last 1-cell of 
1ci((t) if any, the first 1-cell of lc2(o") otherwise, let a be the initial point of 
bo sink(e) or bo sour(e) depending on whether e is the last 1-cell of 1ci((j) 
or the first 1-cell of lc2(o"), and let J be the maximal antichain whose 1-, 
and bounded 2-cells are the cells whose sinks are in the principal filter of the 
sink of e but not their sources. We define inductively a finite sequence of 
pseudotriangulations Gi, G2; ■ • ■ j Gp, and a finite sequence of pseudotriangles 
Ai, A2, . . . , Ap, a G 9Aj, as follows 

(1) Gi = boG(J) and Ai is the pseudotriangle of Gi lying locally to the 
left of the sink bitangent line segment of a if any; otherwise (that is, if 
the sink bitangent line segment of ex is a right-right boundary bitangent 
line segment as in the example of Figure [2T]) we define Ai to be the 
pseudotriangle lying locally to the right of the image under i of the 
source bitangent line segment of a; walking in counterclockwise order 
along the boundary of a pseudotriangle A j , a G (9A j , starting from the 
point a we traverse successively four convex chains Al, A?, A? and A^, 
the last one A^ being the empty chain in case a is a cusp point of Aj; 

(2) If A? is an arc or a boundary bitangent line segment we set i = p 
and we are done; otherwise we define Aj+i to be the pseudotriangle 
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lying locally to the left of the sink bitangent line segment of a in the 
pseudotriangulation Gj+i obtained by flipping clockwise in Gi the first 
bitangent line segment tj of A?. 
Playing the same game with e' the last 1-cell of 1x1(0") if any, the first 1- 
cell of rc2((T) otherwise, we define similarly a sequence of pseudotriangulations 
G'^ , . . . , and a sequence of pseudotriangles A'^ , . . . , A^, . We denote by Lgink 
and Rsink the operators that assign to a bounded 2-cell a the pseudotriangle 
Ap and A^, defined above; note that Rsink and Lsmk are conjugate under l. 
Finally playing the same game with the dual order of the order -< we introduce 
similarly the operators Lgour and Rsour- 

Example 2.10. Figure [22] depicts the sequences of pseudotriangulations Gi, G[ 
together with the pseudotriangles Aj, A'- for the 2-cell a of the visibility com- 
plex of the family of bodies introduced in Example 12.91 : A case where the sink 
of a is a right-right boundary bitangent : Here e is the first 1-cell of lc2(cr) 
and its sink is 623; e' is the last 1-cell of rci((T) and its sink is feg^, p = = 1. 



Example 2.11. Figure [23] depicts the sequences of pseudotriangulations Gi, G'^ 
together with the pseudotriangles Aj, A^ for the 2-cell a of the visibility com- 
plex of the family of bodies introduced in Example 12.81 : here e is the first 
1-cell of lc2 (cr), the sink of e is the bitangent 623, P = 4, the sequence ti, t2> ^3 
of flipped bitangents is the sequence hy^-, bn, byf, e' is the last 1-cell of rci (cr), 
the sink of e' is the bitangent and p' = 1. 

Theorem 2.9. Let a be a bounded 2-cell whose source is an interior vertex. 
Then the pseudotriangles Lsink(o") o.iT'd Rsink(<7) are adjacent along sink(a) and 
their union 'H{a') is a free pseudoquadrangle whose bitangent line segments 
lying on its boundary are the bitangent line segments of Hi.c(o") and Hic(o"), 
possibly augmented with 1 or 2 boundary bitangent line segments in case the 
first or the last or both the first and the last elements of the label of a are 
the infinity symbol. Similarly the pseudotriangles Lsour(c) and Rsour(c) are 
adjacent along sour((7) and their union is T-L{a). □ 

Proof. We denote by the sequence of bitangent line segments of the chain 
A^, j G {1, 2, 3, 4}, i = 1,2, . . . ,p. One can easily check that the sequence Aj 
is well-defined, finite, and that 




Figure 22. 
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Figure 23. 



(1) 



A} 



Rlia), and = H/,(a); 



(2) Af is a prefix factor of A?_,_j^, and Vp = sink{a) ll^^{a) unless the 
sink bitangent of a is a right-right boundary bitangent, in which case 



A2 

p 



A1 is an arc. 



The theorem follows. 



□ 



Theorem 2.10. Let a be bounded 2-cell whose source is a boundary vertex. 
Then the bitangent line segments of the sequence Hrc(o") (resp. Hic((t) ) bound a 
pseudotriangle incident to the boundary of the convex hull along the bitangent 
line segment corresponding to the sink (resp. source) or the source (resp. 
sink) of a depending on whether the source of a is a left-left or a right-right 
bitangent. □ 
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Note that under the assumption that the source of a is a boundary bitangent 
then the pseudotriangle defined by Hj-dcr) is denoted indifferently Rsink(<7) or 
Rsour(f); similarly the pseudotriangle defined by Hic(cj) is denoted indifferently 

Lsink(0") or Lsour(cr)- 

Combining now Theorem 12.61 and the previous analysis we get the following 
key result for our purpose. 

Theorem 2.11. Let a be bounded 2-ceU and let e be a left 1-cell of the 
left boundary chain of a. Then the bitangent ^Phis) leaves the pseudotrian- 
gle Lsijik(cr). (In case the source of a is a left-left boundary bitangent and if e 
is the first left 1-cell of the left boundary of a then (fhie) o,nd the image under 
L of the source of a coincide.) □ 

2.3.2. ABo-pseudotriangles. Let J be a maximal antichain of O free of vertices. 
We use the symbols urj, ulj, drj, dlj for the products urfotopj, ulfotopj, 
drfobotj, dlf obotj where urf(cr) = cif (urf for "upper right face"), ulf{a) = 
a"b, drf((T) = dlf((T) = for a left-left 0-cell or a left 1-cell a, and where 
urf(fj) = ulf{a) = fJi, drf((j) = cJb, dlf (a) = C7f for a right-right 0-cell or a 
right 1-cell a. (Note that urf and dlf on one hand and ulf and drf on the other 
hand are conjugate under ly.) The >l,So-pseudotriangles of a, a G T°{J), are 
denoted £/o{cr) and c^o{cr) and are defined as follows 



(8) ^oia) 
and 
(9) 



(Lsink(dlj(fT)) if CTe7-°|i 2,3}('^) 

1 Rsour(f'") otherwise, 

^o(a) = |^smk(urj(cr)) if o- G -^{i^2,3}i('^) 

1 Lsour(o') otherwise 



where, recall, Fij{J), i,j £ {1,2,3}, denotes the set of o" G J^{J) such that 
botj(cT) G rci((T) and topj{a) £ lcj{a); cf. Theorem 12.31 (For J,K subsets of 
{1,2,3} we set Fjk{J) = Uiejje-ft' ^ote that the operators s/q and 
=^0 are conjugate under l. 

Theorem 2.12. Let J be a maximal antichain of O free of vertices and let 
a G J-°{J) whose source is a left-left (resp. right-right) boundary vertex. Then 
its sink bitangent line segment is a right-right (resp. left-left) bitangent line 
segment and is one of the three sides of £/o{a) (resp. 3^o{a)). □ 

Theorem 2.13. Let J be a maximal antichain of O free of vertices and let 
a G J-°{J) whose source is an interior vertex. Then the source bitangent line 
segment of a belongs to the boundary of £^{){ct), the pseudotriangle lies 
locally on the left side of the source bitangent of a, and the sink bitangent line 
segment of a leaves £/o{a). □ 

Proof. Assume first that a G -7^{2,3}{i,2,3}('^)- Then ^/^{a) = Rsour(<7) and we 
conclude using the definition of Rgour (f) and Theorem 12.101 Assume now that 
fj G J'i{i,2,3}iJ)- Then one can easily check that 

(1) dlj(a) ^ 6; 

(2) bot J (cr) is a left 1-cell of the left boundary of dlj(o"); 
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(3) sour(cr) = (^b, (bot J (cr)) and sink((T) = iph{hotj(a))] 

(4) bosour(a") € H^{dlj{a) ,hotj{a)) (since sour(o") is an interior vertex); 

(5) ^o(t) =Lsink(dlj(fT)); 

from which it follows that the source bitangent line segment of a appears in 
the boundary of £/Q{a). It remains to prove that sink((T) leaves £/o{a); but this 
is exactly the statement of Theorem 12.111 □ 

Thanks to the conjugation relation £/q o t = l o SSq we see that M)((t) 
and 3Sq{(t) are adjacent along the source bitangent line segment of a and 
that the sink bitangent line segment of a joins =2^(0") to ^^^{a). Walking in 
counterclockwise order around the boundary of =2^(0") starting from the tail of 
the source bitangent line segment of a we find successively the convex chains 
=2/0^(0"), s^Q{a), =2^0^(0") (which is reduced to an arc or a boundary bitangent 
line segment), and ^^^{a). Let =2^^^((t) be the concatenation of the chains 
■S^Q^a) and s^Q{a). We know that the sink bitangent line segment of a leaves 
the chain ^^'^{a), that b(botj((T)) is a subarc of an arc, say t~, of the chain 
s^Q^{a) : so we define jaf^'^'^la) to be the suffix subchain of ^/^^{a) starting at 
T~ . Similarly we define the chain ^Q^^(a). By construction the sink bitangent 
line segment of a leaves the chain j2^q^'^((7) and reaches the chain ^q°''((7); 
however one can't use directly these chains to compute the sink bitangent of 
a because neither the source of the first arc of £/Q^^{a) nor the sink of its last 
arc are efficiently computable. 

2.3.3. AB-pseudotriangles. We now assume that the convex bodies of P are 
lifts of convex bodies of a finite family T> of pairwise disjoint convex bodies of 
the topological plane A. The central point of a convex body of T> is defined 
to be the branched point contained in that body if any; otherwise any point 
chosen arbitrarily in its interior. A boundary tangent is a tangent to the convex 
hull of convex bodies of P; all other tangents are said to be interior tangents. 
The backward/forward view of a line is the first /last atom of its label. 

Let s be the initial segment of an interior tangent t to a convex body of 
T> with backward view a convex body 0'. We set P(s) = {0,0'} Q V. Let 7 
be a simple oriented curve in o' U o U s joining the central point of 0' to the 
central point of o with the property that its projection 7' in A is simple, as 
illustrated in the left top and left bottom diagrams of Figure [231 We construct 
a new branched covering A^(]B) of A as follows. We cut C = B U A along 7 
and 7', call the resulting surface and q : Cy ^ C the induced projection; 
we define A^(IB) as the quotient space of by identifying, on one hand, the 
left lift (under q) of 7 and the right lift of 7' and, on the other hand, the right 
lift of 7 and the left lift of 7'. By construction the set Aj,(P) of connected 
components of q~^ (\J{'D U 'D{s))) is a set of pairwise disjoint convex bodies of 
Ag(IB). A similar construction can be done with the terminal segment of the 
interior tangent in place of its initial segment. 

Theorem 2.14. Let s be the initial segment of an interior tangent to a convex 
body o of D with backward view a convex body o' , let a = dlf (s) and let a' = 
dlf(s') where s' is the right lift of s in the visibility complex of \g{V). Then 
sink(a') = sink(a), Hrc(a') = Hrc(a), H+(a',s') = H+((t,s), and H-(a',s') is 
the right-left or right-right bitangent line segment joining o' to o depending on 
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Figure 24. 

whether s is a left or right tangent. A similar result holds for the terminal 
segment of an interior tangent with forward view a convex body. □ 

Proof. By construction modulo the simple observation that the segment s does 
not cross any of the bitangent line segments of Hi.c((t) and Hi(;(cr). □ 

A similar construction is done in the case where the backward view of t 
is the infinity symbol using the left-left boundary bitangent line segment v 
pierced by s and the convex body o' that v leaves, as illustrated in the right 
diagrams of Figure [251 Note that in the case where 4 is a left tangent one can 
have o = o'. 

Theorem 2.15. Let s be the initial segment of an interior tangent whose 
backward view is the infinity symbol, let a = dlf(s) and let a' = dlf(s') where 
s' is the right lift of s in the visibility complex of Xg(T>). Then sink(o"') = 
sink(cr), Rrc{cr') = Urdcr), H+(cr',s') = H+(cr,s), and R~{a' , .s') is the left- 
right bitangent line segment joining o' to o if s is a right tangent, the left-left 
bitangent line segment joining o' to o if s is a left tangent and o ^ d ; the 
empty sequence otherwise. A similar result holds for the terminal segment of 
an interior tangent whose forward view is the infinity symbol. □ 

We turn now to the definition of the >l;S-pseudotriangles. 

Let J be a maximal antichain free of vertices of 0{T>) with a distinguished 
line r(e) G e for each e € £{J) and let s be the initial segment of one of the 
T(e), e G S{J), with the property that the supporting line of s is an interior 
tangent. Let J' be a maximal antichain free of vertices of 0{T>[s)) with a 
distinguished line T(e) € e for each e G £{J') such that J and J' agree along 
the supporting line t of s, that is, the projection of t in A is a T(e') for some 
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e' S £{J')- In that case the hnes (?~^(r(e)), e € £{J) Uf (J'), define a maximal 
antichain of the visibihty complex of Ag(I'), denoted Ag(J) thereafter. 

Let a G Fxy{J) whose source is an interior vertex. We define /i^(B), 
/x+(J), and to be A^(IB), A^(P), A^(J), and dlf(s'), where s is the initial 

segment of r(topj(cr)) and s' is the right lift in /i^(B) of s, if y G {2, 3}; B, P, 
J, and CJ otherwise. We define ^~ to be the conjugate of fi^ under i and we 
set fJL = /i~ o ^+ . Applying twice Theorem 12.141 we get the following theorem 
where we write H^(cr; J), Hf|;(cj; J), H~(cr; J), H[~((t; J), for H^(cr, botj(cr)), 
H[|;(cr,topj(cj)), H7^(cr, bot j(o-)), and ll~^{a,ioY)j{a)). 

Theorem 2.16. Let a € J-xy{J) whose source is an interior vertex, let J' = 
fJ-{J), and let a' = ^{a). Then 

(1) a' G Fxy{J'); 

(2) sink(o-) = sink(cj'); 

(3) H+(c7'; J') = H+(ct; J), H+K; J') = H+(a; J); and 

(4) H~ (cj'; J'), Hj^(cj'; J'), ancf sour(cj') are computable in constant time. 

□ 

As announced in the introduction the ^0-pseudotriangles of cj are the ABq- 
pseudotriangles of /i(cj); they are denoted ^(a) and .^(cj) in the sequel. 

2.3.4. Summary. We reformulate, in preparation to the description of our 
pseudotriangulation algorithm in the next section, Theorem 12.161 using the 
following more adequate notations. For cj € J^°{J) we set 



(10) /Cdi(cj) 



iff^G-^3{i,2,3}(>/); 
{dlj(cj)} iftTG J2{1,2,3}(^); 

{dl J (cj) , ur J o dl J (cj) } otherwise; 



similarly we define /Cur(o') to be the empty set if cj € -^{i.2,3}3(<^); the singleton 
{urj(cj)} if CJ € -^{i.2,3}2('^); the pair {urj(cj), dlj o urj(cj)} otherwise; note 
that /Cur is the conjugate of /Cdi under the shift operator u. Thus Theorem l2.16l 
can be read as follows. 

Theorem 2.17. Assume that the orbits of the sink bitangent line segments of 
the 2-cells of K,^\{a) U /Cur(o') under the operators Pback o.'nd Pforw fl'^e known. 
Then representations of the AB-pseudotriangles ^{a) and =^(cj) are com- 
putable in constant time. □ 

Let G be the digraph whose set of nodes is the set of arcs of the cross- 
section and whose set of arcs is the set of pairs (cj, cj'), cj' G /Cdi(c) U /Cur(<7), 
as illustrated in Figure [25] where we have drawn G on the canonical upward 
drawing of the cross-section. If G is acyclic then, according to Theorem 12.171 
any topological sort of G provides a total order on the set of arcs of the 
cross-section to compute their ^;S-pseudotriangles in constant time per arc. 
Unfortunately the digraph G is not acyclic in general. The following (easy to 
check) theorem will be used in the next section to decide where to break the 
cycles of G in order to be able to compute the ^S-pseudotriangles no longer 
in constant time but in constant amortized time per arc. 
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Figure 25. 

Theorem 2.18. Let a be a left 1-cell of J. Then uri{a) G -^3{i,2,3}('^)) 
uli{a) G •^i{i,2,3}('^)) '^''^d drf((T) = dlf((7) G -^{i,2,3}2('^) U 0. ^4 similar result 
holds for right 1-cell of J using the conjugation relation J^ij o t = l o J^j^. □ 
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3. Our algorithm and its complexity analysis 

We are now ready to describe our pseudotriangulation algorithm. Recall 
that our pseudotriangulation algorithm proceeds in three steps: we first com- 
pute the convex hull of the family of convex bodies, then the cross-section 
of the visibility complex of the family of convex bodies with constraints as- 
signed to a boundary bitangent line segment (not forgetting to add all the 
boundary bitangent line segments computed at the first step to the set of 
constraints), and finally the greedy pseudotriangulation associated with that 
cross-section. The input of our pseudotriangulation algorithm is a finite planar 
family of pairwise disjoint convex bodies together with a distinguished set of 
pairwise interior non-crossing free bitangent line segments of the family: the 
constraints; the family is only given by its chirotope, that is, for all triple of 
indices of the family the position vector of the corresponding triple of convex 
bodies is computable in constant time; equivalently, for any convex body of the 
family the relative counterclockwise circular order of any triple of bitangents 
tangent to that body is computable in constant time; cf. Appendix \^ The 
family is denoted V = {oo,oi, . . . ,On-i}, the underlying topological plane is 
denoted A, and the left-left, left-right, right-left and right-right bitangent line 
segments joining the body Oj to the body o-j are denoted v^, u -t, vi- and vjj, 
respectively. 

3.1. Convex hull algorithm. Let p : B ^> A be a (connected) 4-sheeted 
branched covering of A ramified over the central point of oq, let r be a generator 
of its automorphism group {k, Z4), and let Xi, i G Z4, T{Xi) = Xj+i, be the 
four connected components of the pre-image under IB — )• A of the complement 
in A of a curve 7 C oq U wqi joining the central point of oq to the point at +00 
on the bitangent uqi- The sole lift c of oq is called the central body, and the 
lift of Oj, i 7^ 0, whose interior is entirely included in X^^ A; G Z4, or intersects 
both Xk and X^+i is denoted Oi{k). We denote by T> the family of Oj(A;), 
k G {0, 1, 2}, augmented with the central body. Our algorithm to compute the 
bitangent line segments of the convex hull of T) is based on the following three 
simple observations. 

(1) The set of bitangent line segments of the family T> and the set of pairs 
of crossing bitangent line segments of the family V depend only on the 
chirotope of the family P; 

(2) The convex hull of the family V can be extracted in linear time from 
the convex hull of the family T>; indeed, let S be the counterclockwise 
linear sequence of bitangent line segments that appear in the boundary 
of the convex hull of the family T) starting from the left-left bitangent 
line segment joining the central body to oi(0), one can easily check 
(details are left to the reader) that the first bitangent line segment v 
of S entering an o,j(l) is well-defined, that t{v) appears in the sequence 
S, and that the projection of the factor v . . .w of S where w is the 
bitangent line segment that precedes r(u) in S is the sequence of bitan- 
gent line segments that appear in the boundary of the convex hull of 
v. For example in the configuration depicted in Figure [26] one has S = 

T,ivJ:2Wt{v)T,3 Withf = ^51/, Wt{v) = f4'5''i^5'l", Si = t'0lfl2'y23t'34t^45, 
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S2 = ViiqVq2'V2'3''^3''1' 1 and S3 = t;i//6'W6'2"^^2"3"^3"4"^^4"5"^5"6"^6"0- 

The projection of the factor vT,2'w is the counterclockwise sequence of 
boundary bitangent hne segments; 




X2 



Figure 26. 

(3) The Graham's scan to compute the convex hull of a family of points can 
be generalized to families of pairwise disjoint convex bodies provided 
that the bodies are sorted around a boundary body. Since the central 
body is, by construction, a boundary body of the family D we can apply 
the generalization of Graham's scan we have in mind to compute the 
convex hull of the family V, and apply our second observation to derive 
the convex hull of the family D. 
We now explain our generalization of the Graham's scan. 
We perform a counterclockwise rotational sweep of the 4-sheeted branched 
covering space B with a half-line whose supporting line is a left tangent to 
the central body c at its origin. The lift in sheet X^, A; S Z4, of the half-line 
supporting the bitangent line segment Voi with origin its tangency point upon 
oq is denoted ii{k). The sweep starts at position ii{0) and ends at position 
^1(3). During the sweep we maintain the convex hull of the subset of bodies 
of D that have been entirely of partially swept by the sweeping half-line; to 
this end we keep track of a subset of the bodies that intersect i; therefore 
we update our data structures (to be defined in a second) when the sweeping 
half-line reaches the ii{k): an enter event, and some of the ij{k): a leave event. 
For a given position i of the sweep half-line, we define 

(1) S(^) to be the linear sequence V1V2 ■ ■ -Vk of bitangent line segments 
encountered when walking counterclockwise along the boundary of the 
convex hull of the bodies that has been reached so far by the sweep 
half-line, starting from the left-left bitangent line segment joining the 
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central body to oi(0); the body that the bitangent hne segment Vi 
reaches is denoted (in particular o'^ is the central body); 

(2) v^,(^) to be the bitangent line segment Vji where j' is the minimal 
element of the subset of indexes j (1 < j < A:) such that an infinitesimal 
counterclockwise shift of i pierces the bodies o[ {j <i < k) in the order 
Ofc) o'^-ii ■ ■ ■ I'^'p by construction i pierces either v^,(^), or its successor 
arc or its successor bitangent line segment; 

(3) Q(^) to be the list of bodies o'^o'^_i • • • o'-,. 

For example in the configuration depicted in Figure [27] the set of convex bodies 




X2 X3 



Figure 27. 

entirely or partially swept by the current sweeping half-line I is {0, 1, 2, 3, 4, 5, 6, 1'} 
and one has 

= VQiVi2V22.V2.iVi5V^vVVfiVQQ, 

= 061', 

During the sweep, we maintain the lists Q(^), and the bitangent line 

segment v*(^). Initially i = ^i(O), v*(£) is the left-left bitangent line segment 
^c,oi(o) joining the central body to oi(0), S(^) = bc,oi(o)> V(o),c], and Q(£) = 
[c, oi(0)]. We store Q(£) in a binary search tree. Assume we are to process the 
enter-event £' for the body o successor of the enter or leave event L Let r be 
the rightmost body of Q(^), r^, its predecessor along and r* its successor, 
if any, that is, if r is not the central body. Now, we explain how to update S 
and Q. First we locate o in Q. Assume first that o is to the left of r. Let a and 
a' be its left-hand and right-hand neighbors in Q(^). If o does not intersect 
V = Vaa' then o is included in the convex hull of a and a' and we just ignore 
o: S(f ) = Q(f ) = Q(^) and v*(f ) = v*(£). Otherwise we insert o into 
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Q, and update S: we split S at f , and we regard the two resulting parts as 
stacks of arcs whose respective heads are the arcs contributed by a and a'. 
Then we pop from the left-hand stack until an arc /3, say supported by body 
o', is met such that d is the central body ox v = Vo qi reaches /3. Similarly we 
pop from the right-hand stack until an arc say supported by body o" , is 
met such that v" = Vo"o reaches /3'. Then, we shorten /3 and /?': the source 
of /3 and the sink of (3' are replaced with Vqo' and Vo"o, respectively. Then, to 
build we concatenate what is left of the two stacks, with the arc, say S, 

of do with source Vo"o and sink Vqo' in between. When an arc that follows the 
arc a{i) that v*(£) reaches (included) in is popped, its supporting body 
is removed from Q. If r(i) is removed from Q, and the body supporting the 
predecessor of 6 along T,(i') intersects i' at the right of o, we insert it into Q, 
so that it becomes r{£') instead of o. If r(i) is removed from Q, and the body 
supporting the predecessor of 6 along does not intersect I' at the right of 
o, o becomes r{i'). Assume now that o is to the right of r. We discard o when 
o is included in the current convex hull, that is, if o is included in the convex 
hull of r, and r*. Otherwise, we proceed as in the previous case, except 
that we split through the arc a that v^(£) reaches instead of bitangent 
Va'a (that is, there is one copy of a at the head of both stacks), and a body 
is removed from Q only if an arc it supports is popped from the left-hand 
stack. The body fj, supporting the predecessor of 6 along T,{£') is inserted into 
Q if a(i) has been popped from the right-hand stack and fj, intersects £' at the 
right of o. Finally concerning the leave events only leave events for r need to 
be processed. The processing of those events simply consists in removing r 
from Q . 

Finally we mention that if, instead of working in a 4-sheeted covering, we 
work in a 3-sheeted covering, and lift the bodies accordingly, then the convex 
hull of the lifts still contains the convex hull of the bodies as a factor but 
our criteria to locate efficiently this factor breaks down; see Figure [28] for an 
illustration. 




the plane sheet sheet 1 sheet 2 



Figure 28. Three sheets are not enough. 
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3.2. Cross-section algorithm. As said in the introduction our cross-section 
algorithm is a sweep of the convex hull of the bodies by a half-line whose 
supporting line is a left tangent at the origin of the half-line to a distinguished 
boundary body. The sweep starts at a distinguished boundary bitangent line 
segment leaving the distinguished boundary body (there might be several), 
and during the sweep we construct the cross-section of the visibility complex 
of the family of bodies with constraints assigned to the distinguished boundary 
bitangent line segment. The construction of the cross-section boils down to 
maintaining during the sweep the ordered sequence of forward and backward 
view^ of the lines of free space supported by the sweeping half-line, that is, 
whose second coordinate is the sweeping half-line. However it is asking to 
much of our chirotope because a view might start or (not exclusive) end at 
the endpoint of a constraint, and we have already observed that the relative 
positions of the endpoints of the constraints with respect to the bitangents 
are not completely determined by the chirotope of the convex bodies. To 
overcome this difficulty we are going to embed the views into larger boundary 
curves — called paths in the sequel — that start and end only at the touching 
points of the sweep half-lines with the convex bodies, and not at the endpoints 
of the constraints. So we reduce the problem of computing the cross-section 
to the problem of maintaining during the sweep the ordered sequence of paths 
pierced by the sweeping half-line — each path being represented by a subpath 
of constant complexity, called its window, because contrary to the sum of 
the complexities of the views the sum of the complexities of the paths is not 
necessarily linear. We let the (standard) details of the maintenance of the 
pierced paths by the sweep half-line to the reader and we concentrate on the 
definition of the paths and on the definition of their associated windows. 

Assume without loss of generality that oq is a boundary body and that 
the sweep is done around oq. Let G be the geometric graph union of the 
boundaries of the convex bodies Oj (except oq) and the constraints rj (we 
delete from the set of constraints the constraints incident to the body oq). 
Let nii and Mi be the touching points with the body Oj of the left-left and 
left-right bitangents joining oq to Oj; similarly let Oj and Ai be the touching 
points with the constraint Tj of the left-left and left-right bitangents joining 
Oq to Tj. The points nii and Mi split the boundary of Oj into two arcs gi 
and di where by convention di joins rrii to Mj when walking counterclockwise 
around the boundary of Oj. The arcs gi and di are oriented from nii to Mi, 
and the constraint Tj is oriented from Oj to Ai. This turns the graph G into 
a directed acyclic graph that is monotone with respect to sweeping half-line 
coordinate. We construct a new geometric graph G' in two steps: for every 
body Oi we firstly add to G the chords of Oj joining (1) nii to Mi, (2) nii to 
every aj € doi, and (3) every Aj G doi to the first € doi that follows, if any; 
Aj to Mi otherwise — and secondly we delete the body boundary edges of G. 
(See Figure [29] for an illustration.) The unique maximal monotone path in G' 



^The views are the connected pieces of the boundary of free space cut at its cusp points- 
there is one cusp point per endpoint of constraint — and at tlie contact points of the left-left 
and left-right bitangents joining the distinguished boundary body to the other bodies of the 
family. 
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whose first constraint atom is Tj is denoted p{Tj); similarly the path reduced 
to the chord joining nii to Mj is denoted p{oi). The set of paths p{Tj) and p{oi) 
is denoted V. 




Figure 29. Construction of tlie geometric graph G' . 

We now introduce the window of a path assigned to a sweeping half-line. 
Let p be a path of V and let £ be a sweeping half-line piercing the path p, let 
o' be the first body pierced by i among the bodies — contributing for a chord 
to the path p — lying at the right of the path p and let c' be a chord of the 
path p supported by o'; similarly let o" be the last body pierced by i among 
the bodies lying at the left of the path p and let c" be a chord of the path p 
supported by o". We define a set of at most six chords j = 1,2,3, as 

follows: (1) ei is the highest chord of the path p supported by a body below 
i, 62 is the highest chord of the path p below c' supported by a body pierced 
by i and following (strictly) o', and €3 is the highest chord below c" supported 
by a body pierced by £ and preceding o"; and similarly (2) e'l is the smallest 
chord of the path p supported by a body above i, €2 is the lowest chord above 
c' supported by a body pierced by £ and following (strictly) o', and 63 is the 
lowest chord of the path p above c" supported by a body pierced by i and 
preceding o" . The window of the path p at ^ is denoted f{p,£) and is defined 
as the subpath of p defined as the intersection of the paths p, p+{ej) and 
P-{e'j), j = 1,2,3, where for any chord c of the path p, the subpath p+{c) is 
the suffix subpath of p starting at the constraint following c and similarly the 
subpath p_ (c) is the prefix subpath of p ending at the constraint preceding c. 
A simple case analysis leads to the following theorem. 

Theorem 3.1. Let p be a path ofV pierced by the sweeping half-line i. Then 
the subpath f{p,i) of p is well-defined, contains at most two chords, and is 
pierced by the half-line £. □ 

It remains to recall the definition of the views and to embed the views into 
the paths. Let S be the surface obtained by cutting the free part of the convex 
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hull of the bodies along the constraints — note that S is not necessarily con- 
nected. The cusp points of dS (one cusp point x per endpoint x of constraint) 
and the extreme points mi,Mi of the bodies induce a natural decomposition 
of dS into monotone convex paths Vi, called the views. There are two views 




Figure 30. The views. 

per vertex Oj — a view B{ai) whose first atom is the constraint leaving dj and 
a view A{ai) whose first atom is an arc leaving dj — and also two views per 
vertex rrii — a view F{mi) whose first atom is a subarc of gi leaving rrii and a 
view B[mi) whose first atom is an subarc of di leaving mj, see Figure [501 for 
an illustration. Let Tj be a constraint leaving the convex body Oj. The view 
B[ai) is assigned to the path p(rj) and the view A[ai) is assigned to the path 
Pi'^i) where t[ is the constraint leaving Oj that follows Tj, if any; to the path 
p{oj) otherwise. The forward view F{mi) is assigned to the path pijj) where 
Tj is the first constraint leaving gi if any; to the path p{oi) otherwise. Similarly 
the view B{mi) is assigned to the path p{Tj) where Tj is the first constraint 
leaving di if any; to the path p{oi) otherwise. Note that exactly two views are 
therefore assigned to a path. 

3.3. Greedy pseudotriangulation algorithm. The input of the greedy 
pseudotriangulation algorithm is a cross-section r(J) of the visibility com- 
plex of the set of convex bodies P and the restriction of the sink bitangent 
line segment operator to the subset of 2-cells of that cross-section whose sink 
bitangent line segments are boundary bitangent line segments. 

We denote by <^ the partial order on J-" = J- {J) defined by a ^ cj' if there 
is an edge-path in r( J) with source cr and sink a' (recall that the cross-section 
is acyclic). Let now < be a partial order on T compatible with <C on the sets 
•^{i,2}3 a-iid Fi2-, compatible with the dual order <C* on the sets F^[i^2} and 
F21 , and such that 

(11) J33 < -^3(1,2} U -^{l,2}3 < -^22 < ^12 U T21 < J^ll- 

Our algorithm maintains the directed graph T(J), the restriction sinkj- of the 
sink operator to J and the restrictions to B[J) = i?forw(^) U -Bback(^) of 
the operators Pforw and Pback when J describes a maximal chain of down-sets 
in the interval [0,^] (here i?forw(<7) and i?back(f) are the orbits of the sink 
bitangent line segment of a under Pforw and Pbacki respectively). 

Let i7 be a down-set of (J-", <). Let cr be a minimal element oi J-\J whose 
sink bitangent line segment t is not a boundary bitangent line segment. We 
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explain how to compute B' = B{J' U {c}) from B = B{J). Thanks to the 
conjugation relations J-ijOt = ioTji it is sufficient to examine the cases a G Tij 
with i < j. Assume first that a ^ J-'is. In that case one can easily check that 
B' = Bu{t} and that linked representations ^{a) and ^{a) are computable 
in constant time (cf. Theorem 12.171 and Theorem I2.18|) . So it remains to 
explain how to compute t efficiently. We postpone this point to the next 
paragraph. Assume now that a G J-is. In that case B U {t} might be a proper 
subset of B'; so our goal is now not only to compute t but also its orbit under 
Pback (its orbit under pforw reduces to t since t is right-right of left-right). We 
proceed as follows. Let o"i,(T2, . . . ,crk+i with k > 1 he the sequence of 2-cells 
defined by <ti = ct, o-j+i = dlj((jj), a2,...,ak G -F12, and fjfc+i G -^{2,3}2, or 
= 0. Assume to fix the ideas that cJfe+i G J-22, that its forward view is a body, 
and that the backward view of the cTj is a body (the other cases can be treated 
similarly). In that case dlj((Tfc+i) G -7^{i,2,3}3 belongs to J', the bitangent line 
segments of the sequence H^((Tfc+i; J) are bitangent line segments of B, and 
a linked representation of £/{ak+i) is computable in constant time. Assume 
for a moment that the 2-cells urj((T,) (2 < i < A; -|- 1) are in 7^33. In that case 
the SS{(Ti), i = 2, 3, . . . , -|- 1, are computable in constant time. Therefore one 
can compute successively the sinks of the as the bitangent line segments 
joining the (fij) to the SS{(Ti) for i = A; -|- 1, /c, . . . , 2, 1. Now we drop the 
assumption that the fjj are elements of J-33. Let Sj be the terminal segment 
of a line ti G topj(cJi), let be the right lift of Sj in the visibility complex 
of o • • • o \g^^^[V), and let J be the antichain Xg^ o ■ ■ ■ o Asj.^-^(J). Let 
ai = drf(s9 {i = 2,3, . . . , k + 1) and let di = ulf(s2). Using similar arguments 
to the ones given in the proof of Theorem 12.141 one can prove that the orbit of 
the sink bitangent line segment of ai under pback coincides with the orbit of the 
sink bitangent line segment of ai under Pback! furthermore urj(CTj) G ^33 (J), 
2 < i < k. Therefore one can proceed exactly as in the case urj((Tj) G J-33 but 
with the CTj. The cost of this computation is a big-0 of k plus the size of the 
prefix part of H^(dlj((Tfc+i); J) involved in the computation, which sum up to 
a 0(n) when a ranges over J^i^. 

It remains to explain how to compute the bitangent that leaves £/ (a) and 
enters SS{a) in 0(1) amortized time. We will charge the cost on the bitan- 
gent line segments of the greedy pseudotriangulation G( J) associated with the 
cross-section. Recall (cf. [37]) that a bitangent u G G(J) is said to be right- 
minimal if u is the ^-minimal element of the set of bitangent line segments that 
appear in the boundary of the pseudotriangle of G( J) lying locally to the right 
of u; in that case this pseudotriangle is independent of J and is denoted R(u). 
Similarly a bitangent u G G( J) is said to be left-minimal if u is the ^-minimal 
element of the set of bitangent line segments that appear in the boundary of 
the pseudotriangle of G( J) lying locally to the left of u; in that case this pseu- 
dotriangle is independent of J and is denoted L(n). Clearly the sum of the 
complexities of the R(u) and L(n) for u right or left-minimal in G( J) is linear 
in the size of G(J), that is, a 0{n). Let 6162 . . . and e'^e^ ■ ■ ■ e'^/ {k, k' > 1) 
be the sequences of (oriented) arcs of the chains ^""^^{a) and i2/'<=^(o-), and let 
Cm and e^, be the arcs that t enters and leaves. The source and the sink of 
Ci are denoted Vi and f j+i. Note that for i = 2, 3, . . . , A: the bitangent Vi is an 
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atom of the chain SS'^icf). For i = 2, . . . , A;, let Tj be one of the two pseudotri- 
angles of size two adjacent to ^{a) along the bitangent Vi and let Ui be the 
bitangent joining Tj to =^(ct); it is convenient to denote by ui the source of 
/i(cr). Clearly Uj enters SS^ia^SS'^ia^ — which has constant complexit}!! — and 
Ui is computable in constant time. The bitangent line segments Uj decompose 
the pseudopolygon ^(a) U T2 U Ta U • • • U — which reduces to ^{a) in case 
k = 1 — into k pseudotriangles Ri, . . . , Rk where Ri lies locally to the left of 
the bitangent Ui {1 < i < k). Note that the number of bitangent line segments 
lying in the boundary of Ri is a 0(1) (except for Ri in the case a G J-'u). We 
define inductively the sequences n* (i = 1, . . . , /c*) as follows: 

(1) RI = £/{a) and we define as the bitangent joining RI to Ri, 

(2) if u* enters the arc Cj or the arc then k* = i and we are done; 
otherwise we define R*^i to be the right pseudotriangle of u* in the 
pseudoquadrangle R* U Ri and we define u^^i to be the bitangent 
joining R*^^ to Ri+i. 

Clearly the sequences Ri,u* are well-defined; R* and Ri are adjacent along 
the bitangent Uj; t is the last term of the sequence of u* (that is, t = n^*); u* 
leaves one of the arcs, say e^. , of the sequence e'^ . . . e^,; ml > nt^ > • • • > m' . 

Consequently t is computable in time 0{m + mj). Let now t' be the sink of 
top(cr). A simple case analysis shows that either t = t', m = 1 and m' = m^; 
OT t ^ t', m > 2, t' is a bitangent of =^(cj), t' is left- or right-minimal in G(J) 
depending on whether ^{a) lies locally on the left or on the right side of t' , 
and m and — m' are bounded by the size of L(i') or R(t') depending on 
whether ^(a) lies locally on the left or on the right side of t' . Its follows 
that X](7gJ'°{j){"^('^) + "^i('^)} = 0{n). This prove that t is computable in 
constant amortized time. 



This property is no more valid in the presence of constraints if the family of bodies is 
not well-constrained. 
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4. Conclusion and open problems 

We have presented an efficient algorithm to compute a pseudotriangulation 
of a planar (well-constrained) finite family of pairwise disjoint convex bodies 
with constraints presented by its chirotope. Its running time is a O(nlogn) 
where n is the number of convex bodies of the family. We expect to prove 
its practical efficiency in a forthcoming implementation. The design of our 
pseudotriangulation algorithm relies on an extension of the theory of pseu- 
dotriangulations and visibility complexes to the setting of branched coverings 
of topological planes. We conclude by a sequence of open problems raised by 
this work. Can we extend our pseudotriangulation algorithm to not necessarily 
well-constrained family of convex bodies keeping the same time bound on its 
running time? Can we sweep efficiently the dual arrangement of a finite planar 
family of convex bodies presented by its chirotope using only linear space? Ac- 
cording to the main result of the paper this problem boils down in 0(?7-logn) 
time to sweeping the arrangement of the dual pseudolines of the pseudotrian- 
gles of a pseudotriangulation of the family of convex bodies; this arrangement 
is an arrangement of pseudolines with contact points to which it is tempting to 
apply the topological sweep method of Edelsbrunner and Guibas [2] ; however 
the fact that the chirotope of the pseudolines is not computable in constant 
time (since the pseudolines do not have constant complexities) prevents a naive 
application of the topological sweep method. Can we enumerate the pseudo- 
triangulations of a planar family of convex bodies presented by its chirotope 
in sublinear time per pseudotriangulation using only linear space? Does the 
complex of pseudotriangulations of a planar family of pairwise disjoint convex 
bodies is a (shellable) sphere? a matroid polytope? a polytope? It is known to 
be a polytope for families of bodies lying in the classical real affine plane j41j . 
Finally we observe that these questions can be asked not only for families of 
convex bodies of topological planes but also for families of convex bodies of 
branched coverings of topological planes. 

Acknowledgments. The authors thank the referees for pointing several rele- 
vant references and for their helpful comments and encouragements to improve 
the intuitive presentation. 
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Appendix A. Duality in topological planes 
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Figure 31. A topological plane. 



Topological planes. A topological plane on is a surface homeomorpliic to 
M? endowed with a topological point-line incidence geometry whose line space 
is a subspace of the space of pseudolines of the surface satisfying the so-called 
linear space axiom: any two distinct points p and q are contained in exactly 
one line q, their joining line. It is known that the line space of topological 
plane is a open crosscap, that is, a surface homeomorphic to MP^ with one 
point deleted. A topological plane is called affine if for every point-line pair 
{p, I) there exists a unique line k through the point p that coincides with the 
line / or has no point in common with it; the line k is called the parallel to I 
through p. Classical examples of topological planes are 

(1) the real affine plane or classical topological plane defined as the affine 
topological plane whose set of points is and whose set of lines is the 
set of curves with equations y = ax + b, a,b € and x = c, c € M; 
the real affine plane is the default plane in the field of discrete and 
computational geometry; 

(2) the Moulton planes obtained by starting from the real affine plane and 
replacing the lines of the real affine plane with negative slope by the 
kinked lines {{x, y) | ?/ = oj; -|- 6, x > 0} U {(x, y) \ y = akx + b,x < 0}, 
a € M", & G M, where > 1 is a fixed parameter; 

(3) the Klein model of the real hyperbolic plane defined as the restriction 
of the real affine plane to the open unit disk; and 

(4) the so-called arc planes which include examples of topological planes 
with no embedding into any affine topological plane, e.g., start with 
the real affine plane and replace the lines with positive slope by the 
shifted arcs {{x, y) gM? \ y — n = e^"*"}, m, n G M. 

The projective completion of the real affine plane is the unique desarguesian 
projective topological plane on MP^, a characterization due to Hilbert, 1899. 
The Moulton planes are examples of non desarguesian affine topological planes. 
We refer to [9l HSJ [39l [18] for a more detailed background material on topolog- 
ical planes, and for the precise meaning of the terms used above (topological 
projective planes, projective completion, desarguesian, etc.) but not defined 
explicitly. 
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Duality. We work in an oriented topological plane on M?. The term duality 
refers to the map that assigns to a point of the plane the pencil of lines through 
that point, embedded in the line space of the plane, and more generally to the 
map that assigns to a convex body of the plane its set of tangent lines. The 
basic properties of the duality map for points relevant to our purpose are the 
following 

(1) the dual of a point is a pseudoline, that is, a simple closed curve homo- 
topic to a, hence any, generator of the fundamental group of the line 
space or, equivalently, a non-separating simple closed curve; 

(2) the dual family of a family of points is an arrangement of pseudolines, 
that is, a finite family of pseudolines living in the same open crosscap 
with the property that any two cross in exactly one point; 

(3) the isomorphism class of the dual family of a family of points depends 
only on the chirotope of the family of points and vice- versa; and 

(4) any arrangement of pseudolines is isomorphic to the dual family of a 
planar family of points. 

Figure [32] shows the dual arrangements of the five families of three points real- 
izing the five possible chirotopes on a given indexing set of size three depicted 
in Figure [TJ In this figure the line space is represented by a circular diagram 



with antipodal boundary points identified, one marked point (oo) in the role 
of the deleted point, and one arrow on the boundary of the circular diagram 
to indicate which of the two generators of the fundamental group of the line 
space is the one that fits via duality the orientation of the topological plane. 
According to the third basic property of the duality map of points mentioned 
above the chirotope of a planar family of points can be coded by the isomor- 
phism classes of the subarrangements of size three of its dual arrangement. In 
general the isomorphism class of an arrangement of pseudolines 71 , 72 , • • • , 7n 
can be coded by the poset of its cells ordered by inclusion or, equivalently. 





Figure 32. 
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by its side cycles : Ci , C2 , . . . , where Cj is the circular sequence of signed 
indices of the pseudohnes crossed by the side wheel of a sidecar rolling on 7^, 
according to its orientation, with the convention that an index is signed pos- 
itively or negatively depending on whether the crossed pseudoline is (locally) 
directed towards or away 7^. For example the side cycles of the arrangements 
of Figure [35] labeled A and B are 

1 : 2323 1 : 3322 

2 : 3131 and 2 : 1331 

3 : 1212 3 : 1122, 

respectively, and those of A', B' and B" are obtained from those of A and B 
by suitable permutations of the indices. 

In the companion paper |20^ [21] we show that the above properties of the 
duality map for points extend to the duality map for convex bodies in the 
following terms 

(1) the dual of a convex body is a double pseudoline, that is, a simple 
closed curve homotopic to the double of a, hence any, generator of the 
fundamental group of the line space; 

(2) the dual family of a family of pairwise disjoint convex bodies is an 
arrangement of double pseudohnes, that is, a finite family of double 
pseudohnes living in the same open crosscap, with the property that 
any two meet exactly four times, meet transversely exactly four times, 
and induce a cell structure on the one-point compactification of the 
crosscap; 

(3) the isomorphism class of the dual family of a family of pairwise disjoint 
convex bodies depends only on the chirotope of the family of bodies 
and vice-versa; and 

(4) any arrangement of double pseudohnes is isomorphic to the dual family 
of a planar family of pairwise disjoint convex bodies. 

For example the arrangement of three double pseudohnes depicted in the right 
diagram of Figure [33] is a representative of the isomorphism class of the dual 
arrangement of the family of three convex bodies depicted in the left diagram 
of Figure [33] (the two black vertices are the two tritangents of the family). As 




Figure 33. 
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for arrangements of pseudolines an arrangement of double pseudolines can be 
coded by the poset of its cells ordered by inclusion or by its side cycles which 
are defined similarly, except that there are now two side cycles per index since 
a double pseudoline has, contrary to a pseudoline, two sides : a crosscap side 
and a disk side; thus the side cycle of crosscap type of the arrangement of 
Figure is 



and the one of disk type is 



23322332 
3131T3T3 



21122112 



23233232 
33111133 



12122121. 



Observe that for a simple arrangement the side cycle of disk type and the side 
cycle of crosscap type are the negatives of each other. 

Figure [3l] depicts representatives of the 22 isomorphism classes of non in- 
dexed arrangements of three double pseudolines: Each diagram is labeled 
at its left bottom corner with a symbol to name it (of type where a is 
the 2-sequence of its numbers of 2-cells of size 2 and 3 possibly followed, in 
brackets, with the size of the unbounded 2-cell of the arrangement in the case 
where there are several arrangements with the same 2-sequence; and M* 
are mirror images of one another) and is labeled at its right bottom corner 
with the size of its automorphism group; thus the number (118) of simple chi- 
rotopes of families of three pairwise disjoint convex bodies on a given indexing 
set of size 3 can be computed as the sum 

3! 6 6 6 

^ 1 2 3 

fc>i 

where is the number of arrangements of Figure [Ml with group of automor- 
phisms of order k. 




M*5^(2) 1 M36 

Figure 34. 
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Appendix B. Perturbation scheme 



Our purpose of this section is to give a detailed proof of the following (in- 
tuitively clear) theorem. 

Theorem B.l. Let x be a non-simple chirotope of planar families of pairwise 
disjoint convex bodies. Then there exists a simple chirotope x* j computable in 
constant time, such that x o,nd x* have the same set of pseudotriangulations. 

□ 

Our proof is based on (1) the one-to-one and onto correspondence, induced 
by the duality map, between the class of chirotopes of finite planar families 
of disjoint convex bodies and the class of chirotopes of arrangements of dou- 
ble pseudolines [20]; cf. Appendix [Aj (2) the dual characterization of the 
class of pseudotriangulations of a finite planar family of disjoint convex bodies 
of Pilaud and Pocchiola [33]; and on (3) the classical interpretation of sim- 
ple pseudoline arrangements as non redundant primitive sorting networks [25\ 
section 5.3.4], [Ml page 29], [H]. 

By a pseudotriangulation of a double pseudoline arrangement V living in 
an open crosscap we mean a pseudoline arrangement with contact points £, 
that is, a finite family of pseudolines that intersect pairwise in a finite number 
of points of which exactly one is a transversal intersection point, with the 
property that (J £ covers exactly the topological closure of |J V minus the 
first level of the arrangement V with respect to its noncompact face. Figure [551 
depicts a non simple arrangement of three double pseudolines 1,2 and 3 and 
one of its pseudotriangulation a, 6, c, d (red, blue, green and purple in pdf 
color): the vertices that support a contact point are marked with a disk (yellow 
in pdf color); the non simple vertex supports only one contact point; this 
contact point is a contact point between the pseudolines h and c. 




d 
c 
b 
a 




Figure 35. 



Theorem B.2 (Pilaud and Pocchiola [33]). Let A be a finite planar family 
of pairwise disjoint convex bodies, let V be its dual arrangement, and let <I> 
be the map that assigns to a pseudotriangulation of A the arrangement of 
the dual pseudolines of its pseudotriangles. Then <1> realizes a one-to-one and 
onto correspondence between the set of pseudotriangulations of A and the set 
of pseudotriangulations of S/ . □ 
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Proof. For completeness and clarity we repeat the proof of |33j (that consider 
only convex bodies in general position). That $ is well-defined and one-to- 
one and that its range is included in the set of pseudotriangulations of V was 
explicitly observed in To prove that the range of ^ is exactly the set of 
pseudotriangulations of V we use a simple counting argument involving the 
tangency and winding numbers. Let Ai, A2, . . . , A„ be the convex bodies of 
A and let 71 , 72 • • • , 72n-2 be the envelopes of the pseudolines of a pseudotri- 
angulation of V. Our goal is to prove that the Ti are the boundaries of the 
pseudotriangles of a pseudotriangulation of A. This boils down to proving 
that for any point x avoiding the boundaries of the Aj and 71 

(12) y\(x,%) = l^ ifxGAo\Ur=iA^ 
^—^ ' I otherwise 

1=1 K. 

where w{x, %) is the winding number of x with respect to Ti and where Aq 
is the convex hull of the A.j. Introduce in the picture the tangency number 
t[x, 7i) of a point x with respect to Ti as the number of tangents to Ti through 
X or, to put it differently, as the number of intersection points between the 
dual pseudolines of x and T- Similarly introduce in the picture the tangency 
number t{x,Ai) oi x with respect to Aj. The reader will easily check that 

n 2n-2 

(13) Y.*{x,A,)=t{x,Ao)+Y,t{x,Ti); 

i=l i=l 

that t{x, Ai) = or 2 depending on whether x belongs to the interior or to the 
exterior of Aj; and finally that t{x,T) = 1 + 2w{x,Ti). Equation follows 
easily. □ 

Proof of Theorem \B.1\ Let L be a realization of x as a double pseudoline ar- 
rangement. We define x* as the chirotope of an arrangement of double pseudo- 
lines r* obtained T by a sequence of mutations. Let v he a, non-simple vertex 
of r with degree 2fc, > 3. In the vicinity of v the arrangement V is isomor- 
phic to a pencil of pseudolines ^1 , £2 j ■ ■ ■ ; ^fc of the unit disk as indicated in the 
left diagram of Figure [37] where each pseudoline is oriented from left to right 
and where each pseudoline (except £1 and for which the information is not 
relevant) bears a distinguished side — indicated by a rectangle in the figure — 
which corresponds to the Mobius side of the corresponding double pseudoline. 
We perturb the pencil in the vicinity of its vertex into the unique cyclic 
arrangement l[ with the property that the pseudoline i[ contributes to the 
upper or lower envelope of the arrangement depending on whether its Mobius 
side contains the south or north pole of the unit disk, as illustrated in the 
right diagram of Figure [371 Carrying back this perturbation on F we get an 
arrangement F' whose number of non-simple vertices is one less than the num- 
ber of non-simple vertices of F. This perturbation is clearly compatible with 
the isomorphism relation and stable under taking sub- arrangement, that is, 
for any J C I the restriction to J of the arrangement F' is isomorphic to the 
perturbed version of the restriction of F to J: more formally (F') j = (Fj)'. In 
particular the chirotope of the perturbation of F is the perturbation of the chi- 
rotope of F; this proves that x' is computable in constant time. Repeating this 




s s 



Figure 36. 



perturbation at each non-simple vertex of the arrangement F we get a simple 
arrangement F* whose chirotope x* is computable in constant time. It remains 
to show that x ^iid x' (hence x*) have the same set of pseudotriangulations. 
Let Che a pseudotriangulation of F. We construct a pseudotriangulation C of 
F' which coincides with C except in the vicinity of v. Let a be the permutation 
of {1,2,..., k} that maps the index i on the index j defined by the condition 
that the pseudoline of C (or the first level of F) entering v along ii leaves v 
along £j. Now we interpret the arrangement of lines i'^ as a primitive sorting 
network to carry the keys a{i) from the left endpoints of the £i to the right 
endpoints of the ^^-(j) ([26l page 29]); this yields to a touching or crossing status 
to each vertex of the arrangement i[ depending on whether the corresponding 
comparator is feeded with entries in sorted order or not and, therefore, to a 
pseudotriangulation C of F' which coincides with C except in the vicinity of 
V. For example if, in the example of Figure [371 we take 

_ A 2 3 4 5 6 7\ 
^~V4 123657y 

then we sort (in decreasing order) the array [4, 1, 2, 3, 6, 5, 7] with the sorting 
network 
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corresponding to the arrangement ^[ . The sequence of comparisons and arrays 
obtained during the sorting process is then 

-4 — f-4 4 — f-4 4 — 1-6 6 -|-7 7- 

-1 — ♦-l_|-2-*-2— |-6-i-4— |-7-i-6— f— 6- 

-2 2 —1-1 — 1-6 —1-2 —^7—1-4 -|-5 — *— 5 - 

-3 — 1-6 6— 1-1— |-7-l-2— |-5-l-4— f— 4- 

-6 — 1-3— 1-7 7— l-l-|-5— 1-2-|-3— *— 3- 

-5 — 1-7— 1-3 — 1-5 5 —1-1 — 1-3 —1-2 2 - 

-7 — 1-5 5 -1-3 3 3 —1-1 1 1 - 

(the comparators feeded with a pair of indices in sorted order are drawn 
dashed) which yields to a touching status for the vertices u 12 , f 23 , 'i'34 and 
7^56 of the arrangement l[ and to a crossing status for the other vertices, as 
indicated in the following figure where the touching vertices are marked with 
a little (yellow in pdf color) disk. It remains to show that if C" is a pseu- 



N 




Figure 37. 



dotriangulation of V which coincides with C except in the vicinity of f , then 
C = C" . Let A' be a realization of T' as a planar family of convex bodies and 
let A- be the convex body corresponding to the double pseudoline supported 
by By construction the bitangent line segment hij supported by the line 
Vij intersection of the curves and (.'^^ i < j, is free with respect to the A'- 
if and only if j = f + 1 and the ba^i are pairwise interior noncrossing. This 
prove that C = C" since the complex of pseudotriangulations is strongly flag 
connected. □ 
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Appendix C. Computing the visibility graph of a set of line 

segments 

In this section we show how visibihty graphs of finite planar famihes of 
pairwise interior disjoint hne segments (the so-caUed "polygons with holes" in 
the computational geometry literature [TUl [HUl ESI |3] ) fit into the theory 
of visibility complexes of families of pairwise disjoint convex bodies with con- 
straints. To this end we embed the class of chirotopes of families of points into 
the class of chirotopes of families of convex bodies as follows. Define a thin 
chirotope as a chirotope of a family of convex bodies with the property that its 
restrictions to subfamilies of convex bodies pierced by a line are chirotopes of 
pencils of convex bodies (that is, families of disks with same radius and aligned 
centers of the classical topological plane) , and let ( be the map that assigns to 
a thin chirotope the conjunction of its left-left and right-right components — 
we call this way the maps that assign to a triple of indices the subvectors of 
its associated position vector whose entries are those defined with respect to 
the left-left and right-right bitangents, respectively. It is a simple exercise to 
check that the map ( realizes a constant time computable one-to-one and onto 
correspondence between the class of thin chirotopes and the class of chirotopes 
of families of points. Furthermore for any thin chirotope, say realized by the 
family of convex bodies O = {oj}, the four-to-one map r] that assigns to a 
bitangent line segment joining the bodies Oj to oj the line segment joining the 
points C(oi) to Cioj) has an efficiently computable section ^ mapping a line 
segment on an interior bitangent line segment such that the visibility graph 
of the family of points ((O) with respect to a family of constraints K is the 
image under r] of the visibility graph of the family of convex bodies O with 
respect to the family of constraints £,{K). To see this we introduce a permuta- 
tion q = qiq2 ■ ■ - qn of the pi with the property that there exists a line through 
qi that separates the convex hull of the qj, 1 < j < i, from the convex hull 
of the qji, i < j' < n — such a permutation is computable in O(nlogn) time 
as we shall see in the next section — and we define ^ as the map that assigns 
to the line segment joining C(oi) to C(oi) the bitangent line segment joining 
Oi to Oj whose oriented version with respect to the permutation q is right-left. 
The reader will easily check that the section ^ thus defined satisfies the prop- 
erty mentioned above. See Figure [351 for an illustration. Therefore we get the 
following theorem. 

Theorem C.l. The visibility graph of a planar family of n pairwise interior 
disjoint line segments presented by the chirotope of the endpoints of the line 
segments is computable in 0{size of the visibility graph + nlogn) time and 
linear working space. □ 

We should mention that this mechanism of interpretation of visibility graphs 
of families of line segments by visibility complexes of families of convex bodies 
with constraints (that is, the embedding the map r] and one of its sections 
is implemented in the visibility complex package of the CGAL library [2] (see 
also 131 Remark 5]) and explicitly used to design an enumeration algorithm for 
pointed pseudotriangulations of families of points in general position [8]. 
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Figure 38. This figure suggests how visibility graphs of finite pla- 
nar families of interior disjoint line segments fit into the theory of 
visibility complexes of well-constrained families of convex bodies. 
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Appendix D. Computing a pseudotriangulation of a set of points 

In this section we explain how to compute a pseudotriangulation of a planar 
family of points presented by its chirotope. Let pi,p2) ■ ■ ■ ;Pn be a finite planar 
family of n points, let q = qiq2 ■ ■ ■ qn be a permutation of the pi such that for 
any index i there is line through qi that separates the convex hull of the qj, 
1 ^ J ^ from the convex hull of the qj', i < j' < n, and let e be an orienta- 
tion of the underlying plane. Note that such a permutation is computable in 
0(n log n) time : pick an extreme point in O(nlogn) time using your favorite 
convex hull algorithm and sort the remaining points angularly around it. The 
Graham's scan applied to the sequence q boils down to computing explicitly 
the map u = Ug^e ■ {^2, 93, ■ ■ ■ , qn} {gi, 92, ■ ■ ■ , qn-i} defined inductively by 
the relations 

(14) 

where rj is the first natural number r > such that the triangle spanned by 
ordered triple of points qi,u^{qi-i), and ti^+^(gi_i) is counterclockwise or is 
not defined. For example the following table (where we write i for qi) depicts 
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2 
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4 
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8 
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10 


11 


12 


13 


14 







1 





1 
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1 





1 
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u{i) 
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11 


11 


11 


Si 








2 





1 





1 





2 





3 





2 


v{i) 


1 


2 


1 


4 


4 


6 


6 


8 


6 


10 


4 


12 


4 



the map u and its relative v, obtained by reversing the orientation of the 
plane, associated with the family of 14 points pi of the real affine plane with 
coordinates 

(0, 10), (2, 8.5), (5, 0), (5, 3), (5, 8), (5, 14), (6, 7), (7, 2) 

(7, 5), (7, 11), (8, 9.5), (10, 3), (10, 10), (11, 6) 
and the permutation of the pi obtained by sorting them according to the 
decreasing values of their vertical coordinates, as depicted at the left top corner 
of Figure [391 From the map Uq^^ and its relatives Ug-f,, U-g^e, U-g-e, obtained 
by reversing the orientation of the plane or the permutation of the pi or both, 
one can easily deduce, thanks to [371 Theorem 12, Claim 3], not only the 
convex hull of the pi — under the form of the sequence qn,u{qn),u'^{qn) ■ ■ ■ ,qi 
and its relative obtained by reversing the orientation of the plane — but two 
pseudotriangulations (both pseudotriangulations are pointed in the case where 
the family of points is in general position). Indeed let Ug be the set of line 
segments qiu{qi), 2 < i < n, let Vi be the undirected version of the right-left 
bitangent line segment joining the convex hull of the qj, i + 1 < j < n, to the 
convex hull of the qj, I < j < i and let be the set oi Vi, 1 < i < n — 1. 
Then, as illustrated in Figure [321 

(1) Ug is a spanning tree of the family of pi and U^'^ is a pseudotriangu- 
lation of the family of pi; in the case where the points are in general 





Figure 39. Computing a greedy pseudotriangulation of a set of points. 

position the tree is the primal version of the upper horizon tree of 

Edelsbrunner and Guibas |14j : 
(2) U^g induces a decomposition of the convex huh of the family of pi into 

pseudotriangles and pseudoquadrangles, and its completion by is a 

greedy pseudotriangulation of the family of pi (this was first observed 

in [M]; see also [51 Lemma 11]). 
It is also interesting to mention that if denotes the set of line segments 
qiU^'{qi-i), < r < ri,2 < i < n, then T^"^ is a triangulation of the convex 
hull of the family of pi . 
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